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ABSTRACT 

It is generally believed that gaseous disks around supermassive black hole (SMBH) binaries in centers 
of galaxies can facilitate binary merger and give rise to observational signatures both in electromagnetic 
and gravitational wave domains. We explore general properties of circumbinary disks by reformulating 
standard equations for the viscous disk evolution in terms of the viscous angular momentum flux Fj . 
In steady state Fj is a linear function of the specific angular momentum, which is a generalization of 
(but is not equivalent to) the standard constant M disk solution. If the torque produced by the central 
binary is effective at stopping gas inflow and opening a gap (or cavity) in the disk, then the inner part 
of the circumbinary disk can be approximated as a constant Fj disk. We compute properties of such 
disks in different physical regimes relevant for SMBH binaries and use these results to understand the 
gas-assisted evolution of SMBH pairs starting at separations lO""* — 10~^ pc. We find the following. 
(1) Pile-up of matter at the inner edge of the disk leads to continuous increase of the torque acting on 
the binary and can considerably accelerate its orbital evolution compared to the gravitational wave- 
driven decay. (2) Torque on the binary is determined non-locally and does not in general reflect the 
disk properties in the vicinity of the binary. (3) Binary evolution depends on the past history of the 
disk evolution. (4) Eddington limit can be important in circumbinary disks at late stages of binary 
evolution even if they accrete at sub-Eddington rates. (5) Circumbinary disk self-consistently evolved 
under the action of the binary torque emits more power and has spectrum different from the spectrum 
of a constant M disk — it is steeper {vFi, oc z/^^/^) and extends to shorter wavelength, facilitating 
its detection. Our results can be used for understanding properties of circumbinary disks in other 
astrophysical settings. 

Subject headings: accretion, accretion disks — galaxies: nuclei 



1. INTRODUCTION. 

Physics of gaseous disks around astronomical objects, 
also known as accretion disks, has been one of the most 
important topics in astrophysics since the pioneering 
works of Shakura & Sunyaev (1973), Novikov & Thorne 
(1973), Lynden-Bell & Pringle (1974; hereafter LBP74). 
Historically, a lot of attention has been paid to under- 
standing the properties of disks in which the mass accre- 
tion rate M is constant with radius. Some other vari- 
eties of disks have been studied as well, in particular the 
so-called circumbinary disks — gaseous disks orbiting a 
central binary, which can be a stellar binary or a pair of 
supermassive black holes (SMBHs) in centers of galax- 
ies. The latter type of systems has recently attracted a 
lot of attention since the tidal interaction of the SMBH 
binary with the disk results in the loss of the orbital 
angular momentum of the former and its faster inspiral. 
This may help resolve the so-called "last parsec" problem 
(Yu 2002; Lodato et al. 2009) — staUing of the SMBH 
binaries at separations of 10~^ — 1 pc caused by the inef- 
ficiency of stellar dynamical processes at shrinking their 
orbits. Such circumbinary disks are the prime focus of 
this work. 

Depending on the mass ratio of the binary components 
one can have different modes of tidal coupling of the bi- 
nary with the disk. When the secondary-to-primary mass 
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ratio q = Mg/Mp is very small, the secondary cannot 
perturb the disk significantly and migrates through it in 
the so-called Type I migration regime familiar from the 
studies of protoplanetary disk-planet interaction (Ward 
1997). At higher mass ratios'^ the secondary becomes ca- 
pable of clearing gas from the annulus around its orbit, 
switching its orbital evolution into the so-called Type 
II migration regime. As q gets closer to unity the gap 
around the orbit of the secondary turns into a central 
cavity inside of which the binary resides. Numerical sim- 
ulations (MacFadyen & Milosavljevic 2008; Cuadra et al. 
2009) typically show that for q ^ 1 the size of the cavity 
is about twice the semi-major axis of the binary orbit. 

Deposition of the angular momentum of the density 
waves excited by the binary causes non-trivial evolution 
of the circumbinary disk. The seminal study of LBP74 
followed by works of e.g. Lightman & Eardley (1974), 
Lin & Papaloizou (1996) have shown that the behavior 
of disks evolving under the action of both the internal 
viscosity and the external torque can be described by a 
simple equation 
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where A is the external torque per unit mass of the disk. 
Here S and v are the surface density and kinematic vis- 
cosity of the disk, t and r are time and radius, and 

The exact transition between the Type I and II regimes de- 
pends not only on q but also on the viscosity and other disk pa- 
rameters, see Lin &: Papaloizou (1986), Rafikov (2002). 
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I — ri(r)r^ is the specific angular momentum for circular 
orbit. 

In this work, following LBP74 and Lyubarskij & 
Shakura (1987) we reformulate equation ([T]) in a more 
convenient form which allows a very straightforward in- 
terpretation of the steady state disk structure and pro- 
vides a transparent way of understanding the evolution of 
circumbinary disks. This allows us to describe the struc- 
ture and electromagnetic signatures of the circumbinary 
disks around SMBH binaries, their variation in time, and 
the back-reaction onto the orbital evolution of the central 
binary, leading to its inspiral. 

Our work is structured as follows. In ^ we provide 
a general description of the coupled disk-binary evolu- 
tion, derive governing equations, and obtain steady state 
solutions that generalize previously known results. In 
H'S 3.1l we provide the description of the properties of disks 
around SMBH binaries parametrizing them via the vis- 
cous angular momentum flux rather than the mass accre- 
tion rate M . This parametrization allows simple descrip- 
tion of the properties of steady state circumbinary disks. 
In §4] we discuss evolution of circumbinary disks, in par- 
ticular deriving the self-similar solutions to an evolving 
disk structure. All these results are then used in [J5] to 
describe the coupled evolution of a SMBH binary and 
a circumbinary disk, including the self-consistent time 
variation of the disk properties, the orbital inspiral of the 
binary components, and the electromagnetic manifesta- 
tions of the system. Our main results are summarized in 

m 

2. PROBLEM SETUP AND EVOLUTION EQUATION. 

We consider a binary consisting of two point masses 
Mp and Mg < Mp (total mass is Mc ^ Mp + Mg), mov- 
ing around common barycenter on circular orbits (for 
simplicity we neglect the possibility of non-zero eccen- 
tricity). Binary is surrounded by a proghrade, coplanar 
disk (cf. Nixon et al. 2011) which extends to much larger 
distances than the binary semi-major axis ri,. We as- 
sume the eccentricity of the binary to be negligible, even 
though simulations show the possibility of eccentricity 
growth due to the tidal binary-disk coupling (Roedig et 
ah 2011). 

The disk is truncated by the binary torque at the radius 
Tin] the width of the gas-depleted annulus between the 
orbit of the secondary and ri„ is A = | ri„ ~ ri,\ < r?, for 
g <C 1 while A ~ r?, for g ^ 1 (MacFadyen & Milosavl- 
jevic 2008). For simplicity we assume that neither the 
primary nor the secondary have their own disks (this sim- 
plification can be easily relaxed) . The detailed conditions 
for gap opening by a massive perturber in a disk can be 
found elsewhere (Lin & Papaloizou 1986, Rafikov 2002). 
Here we just assume that the binary torque clears out a 
clean central cavity in the disk. In practice this generally 
sets a lower limit on the mass ratio of the binary compo- 
nents q (determined by the local conditions in the disk), 
below which the tidal torque of the secondary is too weak 
to prevent viscous refilling of the gap (or cavity). We de- 
termine the conditions under which such an overflow of 
the disk across the orbit of the secondary (Kocsis et al. 
2012a,b) is possible in ^5 5.2 5.2.61 

External torque A due to tides raised by the central bi- 
nary is expected to be concentrated right at the edge of 
the gap or cavity around the orbit of the secondary. This 



expectation is borne out in calculations of the torque den- 
sity distribution both in uniform disks, where it is gen- 
erally found that A oc \r — rt,\~'^ (Goldreich & Tremaine 
1980; Armitage & Natarajan 2002), and in nonuniform 
disks, where it has been shown by Petrovich & Rafikov 
(2012) that A decays exponentially near the disk edge. As 
a result, already at small separations from r^^ the term 
proportional to A in equation ([1]) can be neglected; this 
equation then reduces to its classical form first obtained 
by LBP74. Effect of the binary torque is then incorpo- 
rated in the solution of this simplified equation via the 
boundary condition (discussed in §2 2.1|) imposed at the 
inner edge of the disk. 

Throughout this work we will assume the potential in 
which the disk orbits to be faithfully represented by a 
Newtonian potential produced by a combined mass Mc- 

Evolution of an accretion disk is best illustrated if we 
characterize the disk at each radius not by its surface 
density E but by the viscous angular momentum flux Fj 
defined as 



Fj = -2nvT,r^ 



dr 



(2) 



where the last equality is for Keplerian disks with 
fl — {GMc/r^y^'^ . This quantity represents the viscous 
torque exerted by the part of the disk interior to a given 
radius r on the external part of the disk, and is thus 
equal to the amount of angular momentum crossing the 
disk circumference 27rr per unit time due to the action 
of viscosity. 

In the absence of external torques mass accretion rate 
through the disk can be directly expressed (LBP74) 
through the divergence of the viscous angular momen- 
tum flux as (note that we take Al > for mass inflow 
towards the center of the system) 

dFj 
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(3) 



motivating us to change the independent variable from 
r to L In this case the evolution equation ([ij with A = 
takes on a particularly simple form (LBP74; Filipov 
1984; Lyubarskij & Shakura 1987): 
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where the function similar to the diffusion coefficient 

o dQ dl 
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is in general a function of both I and Fj because of the 
possible dependence of v oxi Fj. 

Despite the mathematical simplicity of the evolution 
equation Q, which was first recognized by LBP74, over 
the years it has become conventional to study disk evolu- 
tion using the more complicated equation ([1]). We show 
next that use of equation (|4]) allows certain advantages 
over the standard approach, in particular for obtaining 
the steady state solutions for the disk structure. 

2.1. Boundary conditions. 
Using equation ([T]) and continuity equation one finds 
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As we mentioned before, A is expected to be significant 
only in a narrow annulus at the inner edge of the disk. 
We may then assume for simplicity that A = outside 
of some radius ta, which is not too different from (to 
be specific, one can e.g. take ta to be the radius inte- 
rior to which the binary exerts 90% of its torque on the 
disk) . Multiplying equation ^ hy dl/ dr and integrating 
between r-f, and ta one gets 



Fj{rA) = 27r / rAEdr 



Mrirdr, 



(7) 



where we set Fj{ri,) = because S(rb) = by our as- 
sumption of a clean gap. 

The first integral on the right hand side of this expres- 
sion is the total torque that the binary exerts on the disk. 
As long as the orbital evolution of the binary is driven 
predominantly by the tidal coupling to the disk (and not 
due to the gravitational wave emission) conservation of 
the angular momentum ensures that this term is equal 
to 



dLf, Lf, 
2rb 



dt 



(8) 



where Lb = Mc{GMcrbY/'^q/ {I + qY is the orbital an- 
gular momentum of the binary and Vb = drb/dt is its 
inspiral speed. 

In the second integral in equation (jT)) one can write 
M = 2TTY,rvr, where Vr is the radial velocity of the gas, 
and approximate Vr ~ Vb in the annulus between rb and 
rA, thus assuming that gas in this annulus closely follows 
the shrinkage of the binary orbit. Then one can easily 
see that as long as the mass of the secondary Mg satisfies 
the condition 
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^{rb)rb 



Y.nr'^dr, 



(9) 



the second term in the right hand side of equation ([7]) 
is smaller than the first one. This condition is of- 
ten replaced by demanding that the "local disk mass" 
Md = Sr^ in the vicinity of the binary (i.e. at r ~ Vb) 
be less than the mass of the secondary (SC95, hereafter 
SC95; Haiman et al. 2009). 

Whenever we are in the limit Q the viscous angular 
momentum flux in the inner region of the disk can be 
directly related to the orbital evolution of the binary: 



dLb 
~dt 



(10) 



Here we replaced ta with ri„ since the two radii are very 
similar (and also very close to r^). On the contrary, if 
the condition © is not fulfilled the inner part of the 
disk absorbs most of the angular momentum brought in 
by viscous torques and equation (|10|) becomes invalid. In 
this case a more general boundary condition in the form 
(O must be employed. 

Using definition ([2]) and expression for Lb one can 
rewrite equation (jlOp as a formula for the orbital evo- 
lution time of the binary t^v — \d\a.rbl dt\~^: 



Lb 



M, 



2Fj(r„) 67r{l + q)Md' 



(11) 



where we defined the viscous timescale t^, = r'^/v, evalu- 
ated at Vb- Similar expressions (up to a constant factor) 
have previously been quoted in the literature on Type II 
migration in the limit of Md <C Mg, see e.g. Lodato et 
al. (2009), Baruteau & Masset (2012). A subtle point in 
this expression is that both ti, and Md vary in time as 
the disk evolves even if the orbit of the secondary does 
not change appreciably. Calculation of this evolution is 
one of the goals of our present work. 

Previously SC95 have used equation PH)) coupled with 
the condition — as r — >■ r^^ as a boundary condi- 
tion for the problem of the circumbinary disk evolution. 
They have also assumed that a steady state solution for 
the disk structure can always be obtained even far from 
the binary, despite the fact that the radius rb (at which 
the inner boundary condition is imposed) evolves as the 
secondary migrates. Because of this inconsistency their 
solution explicitly depends on time via the dependence 
on rb{t) (and the torque on the binary is found to depend 
on the mass of the secondary) and should evolve on the 
migration timescale of the secondary. However, the latter 
is much shorter than the viscous timescale of the disk far 
from the binary, where this solution is still assumed to be 
valid (by construction, the migration timescale of the sec- 
ondary is of order the viscous timescale at ri„ ~ rb). As 
a result, viscous transport in the disk is unable to com- 
municate information about the changing inner bound- 
ary condition outside the immediate vicinity of the inner 
edge of the disk, implying internal inconsistency of the 
quasi-steady solution derived in SC95. This issue has 
been first noted by Ivanov et al. (1999; hereafter IPP); 
nevertheless the SC95 solution is still being used in stud- 
ies of the disk-assisted SMBH binary evolution (Haiman 
et al. 2009; Kocsis et al. 2011; Yunes et al. 2011). 

In our work we take a different approach and spec- 
ify the inner boundary condition in the form of the 
constraint on the mass accretion rate M. Previously, 
IPP have used boundary conditions similar to ours, even 
though they have formulated them via the asymptotic 
behavior of the disk surface density E. For that reason 
many of our results coincide with their findings. 

At large separations we normally take the disk to be a 
standard constant M disk with the mass supply rate Moo • 
But close to the binary, as r — > 7'i„ , M does not in general 
have to be equal to Moo ■ If the tidal interaction with the 
binary presents a strong barrier for the gas inflow, then 
M(rin) is vanishingly small (Liu & Shapiro 2010). But in 
general tidal torques do not have to completely stop the 
mass inflow: in some situations (e.g. if the gap cleared 
out by the secondary is not deep/broad enough to present 
a serious obstacle to the gas inflow) a fraction of mass 
arriving at the inner edge of the disk can cross the orbit 
of the secondary and be accreted by one of the binary 
components (Kocsis et al. 2012a, b). Alternatively, gas 
can be removed from the circumbinary disk in the form 
of a wind, see ' i'i 3.1 3.1"!^ 15 5.2 5.2.31 To account for this 
possibility we generally use an inner boundary condition 
in the form 



dFj 



dl 



M{r,, 



(12) 



see equation ([3]). Here x 1 is assumed to be constant, 
allowing M{lin) to be less than A'l oo- In practice the 
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value of M{rin) and x is set by the strength of the tidal 
barrier (Liu & Shapiro 2010) and may vary in time as 
the system evolves. 

Clearly, x = implies no gas inflow across the orbit 
of the secondary, and this is the situation that we will 
often consider in this work. For x = 1 secondary does 
not present any barrier to the mass inflow and the disk 
structure reduces to that of a disk with constant M ~ 

Moo- 

2.2. Steady state solution. 

Equation ^ clearly admits a simple steady state so- 
lution (LBP74) 

Fjil)=Fj^o + FjJ, (13) 

where Fjq and Fj i are constants, and / is the specific 
angular momentum. This solution is completely inde- 
pendent of the detailed physics that determines the disk 
properties, since for any (even highly nonlinear) depen- 
dence of Dj on Fj and I the solution (fT3|) still satisfies 
equation ^ in steady state. 

According to equation ([3]) this solution implies M = 
Fj^i also being constant, which seems to suggest that 
equation ((T^ corresponds to the conventional accretion 
disk with constant M. This however is not true is gen- 
eral. 

First description of a constant Af accretion disk was 
provided in a seminal work of Shakura & Sunyaev (1973), 
who explored properties of disks affected only by the in- 
ternal viscous stresses all the way to the central object. 
They have shown in particular that in such disks the lo- 
cal surface density T,(r) at each radius is related to the 
mass accretion rate via 

M = SiTiyj:, (14) 

while the energy loss per unit surface area of the disk 
scales as 

aT^%{r) = ^Mn'. (15) 

Over the years the concept of constant M disks has 
evolved to essentially imply disk properties given by 
equations (ITil) and P3|) . In the rest of this work we 
will call such disks standard constant M disks. 

Let us now consider a steady state solution in the form 
(fT3|) with Fjfl 7^ and Fj^i — M — const. Combining 
equations ^ and ([T^ one finds 

M = 37rz/E - Fj,or\ (16) 
which reduces to only when Fj q — and Fj — 

Ml, and results in a quite different expression for M 
otherwise. In particular, it follows directly from ^ that 
it is possible to have a steady disk with M = as long 
as Fj is independent of radius. In such constant Fj disk 
surface density is related to Fj via 

Fj = Fj,o = 3^i^S(r)/, (17) 

which replaces equation (ITil) . We explore properties of 
such disks in fj3l 

Similarly, the viscous energy dissipation rate in the 
disk per unit radius (and per unit time) dE/dr is given 
by 



Since dE/dr = 4TTraT^ff{r) we can write 



dE_ 

dr 



-Fj 



dn 

dr 
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3 FjQ. 



(19) 



which in steady state described by the solution p3 
yields 



Mn^{r) 



Fj,on 



(20) 



Again, this expression reduces to the conventional result 
only when F/.o = 0. In the case of Af = disk with 
Fj = Fjfi one finds 



Tcff(r) 



/ 3 Fj^/GlTc 
V 8^ a 



1/4 



-7/8 



(21) 



so that Tcg{r) increases towards small radii more steeply 
than in a standard constant M disk (for which Tcs{r) oc 
^-3/4y This result was first obtained by SC95. 

Equation ([T51) predicts non-zero T^g for a disk with 
Fj = const, even though M — in such a disk. Since 
there is no inward mass flow in this disk and correspond- 
ing release of potential energy is absent one may nat- 
urally wonder where does the energy emitted from the 
disk surface come from. Integrating equation (jl8p be- 
tween arbitrarily chosen inner and outer radii and Tq 
one finds the global energy release between these radii 



Eir, < r < r„) = Fj [17(r,) - fliro)] , 



(22) 



(18) 



where we have used the fact that Fj is independent of 
r. Thus, the global rate of energy generation by viscous 
dissipation is equal to the work done on the disk by the 
viscous stress at its inner and outer edges. The latter 
ultimately provides the energy source for the radiation 
from the disk surface. 

To summarize, the assumption of a steady state does 
not necessarily require M and Tcff in a disk to be given 
by standard equations and (|15p — these equations 
are sufficient but not necessary characteristics of steady 
disks. The most general steady state is in fact described 
by equation and then it follows that M and Teff must 
be given by equations p^ and ([20l) . In the rest of the 
paper we explore the properties of such "non-standard" 
disks. 

3. DISK PROPERTIES AS A FUNCTION OF Fj . 

It is conventional to describe the structure of steady 
state accretion disks using the mass accretion rate M, 
assumed to be constant with radius, as a free parameter 
(e.g. Shakura & Sunyaev 1973). Such calculations uni- 
versally assume that relations ([14]) and psp hold true, i.e. 
that Fj^Q = in equation p3|) . Having shown in ^2 2.21 
that in general steady state does not require Fj^q = we 
now revise these scalings by expressing disk properties in 
terms of Fj , which may in general be a function of both 
time and radius r or specific angular momentum I, rather 
than M. When an additional assumption Fj — const is 
made these relations naturally describe the properties of 
M = disks, see 

In our calculations we will always make an assumption 
of an optically thick disk, but it is easy to extend them to 
optically thin disks as well. Because of the large range of 



radii spanned by disks in some systems, the conditions 
in them are expected to vary with radius, resulting in 
transitions in the opacity behavior. To account for this 
in Appendix |X] we derive a set of scahng relations ap- 
plicable to the gas pressure dominated disks with power 
law opacity behavior, and apply them in §3 3.11 to disks 
around SMBH binaries. 

3.1. Circumbinary disks around SMBH binaries. 

Structure of the inner region of a circumbinary disk 
around a SMBH binary may be strongly affected by the 
radiation pressure. We parametrize its role by a di- 
mensionless ratio of the gas pressure to total pressure 
l3 = pg/p = pg/{pr + Pg), where Pg = pkT/n is the gas 
pressure and Pr — is the radiation pressure. In 

the radiation pressure dominated case /? <C 1, while in 
the gas pressure dominated case 1 — /3 <C 1. 

Behavior of viscosity in the radiation pressure domi- 
nated fluid is not well understood at the moment. There 
is still an ongoing debate whether it is determined by the 
full pressure or just the gas pressure in the disk. For that 
reason we use a prescription motivated by the conven- 
tional a-parametrization of Shakura & Sunyaev (1973) 
and accounting for the two possibilities in a convenient 
form (Goodman 2003): 



(23) 



Here b = corresponds to kinematic viscosity propor- 
tional to the full pressure p, while b = 1 corresponds to 
the case when only the gas pressure Pg determines the 
viscosity. Lightnian & Eardley (1974) suggested that ra- 
diation pressure dominated disks with b — are ther- 
mally unstable but recent numerical work by Hirose et 
al. (2009a, b) does not show this to be the case. 

Disks around SMBH binaries are heated by internal 
viscous dissipation and also by energy deposition of the 
density waves launched by the central binary. For sim- 
plicity we do not consider here the latter contribution 
(its role has been studied by Lodato et al. 2009) and 
assume that disk heating is described by equations (|18p 
and (fT9|) . Using equation (|211) we estimate the effective 
temperature of the disk as 



Teff(r): 



3 FjVGMc 



1/4 



-7/8 



(24) 



a.i X 103K ^^]/:a</V,^/«, 



where Me,„ = A4/(10"Mq), r„ = r/(10"pc), and Fj^^ = 
-F'j/(10"erg). This estimate assumes a particular value 
of Fj, which we motivate in ^ ^3 3.1 3.1.41 

In our treatment of the vertical radiation transfer in 
the disk we follow Goodman (2003) and relate the mid- 
plane disk temperature T to TcS via 



rpl 



-T, 



off J 



(25) 



typical for optically thick (r ^ 1) disks, where r 
is the optical depth. Also, 

S = 2ph, (26) 

where p is the characteristic disk density, h — Cg/rt is the 
disk scale height and = {vIp)^^'^ is the sound speed. 



determined by the total pressure p. In the innermost 
region of the disk where the transition between the ra- 
diation and gas pressure dominated regimes occurs we 
assume opacity to be dominated by electron scattering, 
i.e. K = Kes ~ 0.4 cm^ g~^. Further out, in the gas 
pressure dominated part of the disk, k is determined by 
the free-free opacity. We separately consider all these 
regimes below. 

Combining equations (O, (HH), ([25]), & ^ one 

finds the following equation determining the value of /3 
for a disk with opacity dominated by electron scattering: 

- - iGM^)"\..n (27) 



= 2V 



aFj 



29/2 



(1 - /3)i" 

From that one can find the distance r"'^/^'''' of the tran- 
sition between the radiation and gas pressure dominated 
regimes by setting /3 = 1/2 in equation ([27]): 

-. 2/29 



rad/gas _ 



^3.8 X 10-^ 



(28) 



pc 



-| 2/29 



M, 



1/2 



where Q!„ = a/ 10", /i„ = nrup, i.e. ^ in this formula is 
normalized by the molecular weight of fully ionized H. 
Note the extremely weak dependence of r'^^'^/s^^ on the 
binary mass M^. 

3.1.1. Radiation pressure dominated regime. 

Interior to r'^^'^/s'^^ disk is radiation pressure dominated 
and by assumption k = k^s. In this case one finds 

h{r) _ 1 Kes Fj 



r 2n c(l - /3) nr^ 



(29) 



a.l X 10-3 p^^^^^jJ/2^J^/2^ 



irrespective of the value of b in equation ([25)) . At the 
same time the scalings of S(r) and T(r) explicitly depend 
on the viscosity behavior: for 6 = 

Air 

S(r): 



T(r) = 



3 ukI^Fj 
2.1 X lO^g cm-2 

(GM,)V2c2ni/4 



(y■_^ F J T^r^r 



J,50' -2' 



-3/8 



(30) 
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while for 6 = 1 



E(r) 



= 1 




r 24 


(-) 


3^5 7r3 


\ak) 



aF^j 



1/5 



-9/10 



(32) 



■ 6 X 10 g cm 



1/5 



r(r) = 



1 /i(GM,)i/2^, 



-.1/5 



247r2 k 



-11/10 



(33) 



■1.2 X lO^K 



Mo. 5 



Ml','Fl,, 
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Note that a radiation pressure dominated disk with 
Fj = const (and M = 0) has h{r) cx r"-'^/^ (as op- 
posed to the case of a conventional M — const disk for 
which h{r) — const), i.e. the disk puffs up as r decreases. 
Surface density shows dramatically different behavior de- 
pending on the value of b: it rises with r when viscosity 
is proportional to the total pressure (5 = 0), but drops 
with r when v scales with the gas pressure. Finally, T{r) 



is independent of Fj when b - 
5=1. 



but scales as oc F^^ for 



3.1.2. Gas pressure dominated regime with k = k^s- 
Outside of r'^^'^/s^s pressure dominates, meaning 
that /? — 1 and v — ac^/fl, but initially opacity is still 
determined by electron scattering, so that k = k^s. It is 
easy to see that in this case the midplane temperature 
and surface density runs in the disk should be the same 
as in the radiation pressure dominated case with 6=1: 
derivation of T(r) and S(r) involves only the equation 
of vertical radiation transfer ([25]) and the definition ([2]), 
which are identical in two cases because the viscosity and 
opacity are the same when 5=1. Thus, the behavior of 
E(r) and T(r) in the gas pressure dominated case with 
K = Kes is given by equations (|32|) and (p3| correspond- 
ingly. 

However, the scaling of the aspect ratio of the disk with 
r is different from h/r given by equation (|29[) since unlike 
the case studied in §3 3.1 3.1. li the vertical support is now 
provided by the gas pressure. As a result one finds 



h{r) 



1 fk 



l^esFj 



167r2 a<j{GMcfl'^ 
6.6 X 10"^ 



-1 1/10 



1/10 



.-1/20 



-1/20 



(34) 



Note that h/r goes down with r meaning that the disk 
is not flared. However, this dependence on r is so weak 
that the aspect ratio is essentially constant with radius. 

3.1.3. Gas pressure dominated regime with n = Kff. 

At even larger distances opacity in the disk is domi- 
nated by the free-free opacity Kff = Hf f^pT^'^^'^ , with 
Kff^o = 8 X 10^^ cm"i g"2 K'^/^. Using equations ((32|) 
applicable for the situation described in t ^3 3.1 3.0] and 
the behavior of the midplane temperature inferred from 
equation (j34p one finds the transition between the two 
opacity regimes to take place at 
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Outside this radius parameters of the disk scale as 



S(r): 
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, (37) 
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see equations (|A2I) - (|A4I) in Appendix 1X1 In this opacity 
regime disk is only weakly flared and hjr is essentially 
constant with radius. 

3.1.4. Gharact eristic values of Fj. 

We now try to motivate the characteristic values of the 
viscous angular momentum flux in disks around SMBH 
binaries using different arguments. In the case of con- 
stant M disks a standard way of providing the character- 
istic value of M is through the Eddington mass accretion 
rate 

AnGMr^ „ „ _ -iMc,7 



MEdd = 



0.2 Mq yr~ 



£0.1 



(39) 



where e = O.leo.i is the radiative efficiency. 

In the case of constant Fj (or M = 0) disk one can 
also formulate the Eddington limit on the value of Fj, 
which is reached when h/r ~ 1. When this happens ra- 
diation pressure puffs up the disk to such extent that 
it becomes geometrically thick. Further increase of Fj 
should result in mass loss caused by the radiation pres- 
sure driven wind. Another possibility for h/r 1 is the 
mass overflow across the orbit of the secondary which 
may become possible since high above the midplane the 
potential of the secondary is weaker, potentially allowing 
the gas to cross the orbit of the secondary high above the 
disk midplane''. Both these processes result in mass loss 
for the circumbinary disk at its inner edge, which can 
have important implications for the orbital evolution of 
the central binary as we demonstrate in 21 

Unlike the situation with A/Edd in conventional con- 
stant M disks, the Eddington limit on Fj found from the 
condition h/r = 1 depends on the value of r at which it 
is evaluated, see equation ([29]) : 



F,iEdd{r) = 2tt 



(GM,) 



1/2. 



.3/2 



(40) 



It is clear that this condition is most constraining at the 
inner edge of the disk ri„ ^ rf,, where is the semi-major 
axis of the SMBH binary. 

Obviously, one needs to invoke additional considera- 
tions to pick a particular value of (or rf,) at which 
^j,Edd is to be evaluated. Here we assume that this 
critical value of r.m is such that at rb ~ rm the orbital 
decay timescale of the binary due to emission of gravita- 
tional waves tew is equal to the characteristic timescale 
on which the orbit of the binary shrinks due to the tidal 

* For g <g; 1 puffing up of the disk may cause the overflow of the 
secondary orbit even when h/r < 1. 
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coupling to the circumbinary disk tj. The logic behind 
choosing this condition is that the disk then stays sub- 
Eddington all the way until the point when the GW emis- 
sion becomes more important for the orbital evolution of 
the binary than its tidal coupling to the disk. Beyond 
this point the disk becomes super-Eddington and starts 
losing mass in a radiation pressure driven wind, reducing 
the disk torque. However this does not affect the binary 
inspiral since the orbital evolution is no longer sensitive 
to the tidal torque. 
Equating 

^qs c \RsJ 

(here qs = 4q/{l + and Rs = 2GM^/c^ is the 
Schwarzschild radius of the black hole with the mass 
equal to the total mass of the binary Ale) and 



lO^s 10-=^ 10-" 10-3 0.01 0.1 



GW 



in) 



(41) 



t.j[n) 



d In 7'f 



dt 



Lb 
2Fj 



qs McVGM~Fb 



(42) 



(here Lb = ((js/4)Mc(G'Mcrf,)i/2 is the total orbital an- 
gular momentum of the binary) one finds 



^25/25 Fj J ■ 

Plugging this value of r.m into equation pO|) one finds 
that the Eddington value of the viscous angular mo- 
mentum flux in M = disk based on the condition 
towinn) = tjinn) is given by 



F 



J,Edd : 



2103V (GAf,)i4M3Q6ii/io 



'10 



51 



53 C^K^j 
,,17/10 3/5 

erg / qg' . 



(44) 



^c,7 'JS ■ (45) 
This argument justifies the adoption of a characteristic 
value of Fj — 10^" erg (corresponds to Mc = 10^ Mq, 
q — 0.005) in our numerical estimates. 

One can come up with other ways of choosing critical 
Tin or the characteristic value of Fj. In particular, one 
may demand the disk to stay sub-Eddington until the 
point when ^gw becomes equal to the viscous timescale 
at rim after which the binary orbit shrinks faster than 
the viscosity can refill the central cavity (Milosavljevic 
& Phinney 2005). This happens at considerably smaller 
value of Tin than that given by (|43|) and implies lower Fj . 
However, (1) at this point tidal torque is already com- 
pletely negligible compared to the angular momentum 
loss due to the GW radiation, and (2) such condition 
would depend on the poorly understood value of b in 
equation (|23p. Thus, we avoid this way of constructing 
Eddington limit- based estimate of Fj. 

One can also evaluate Fj based on arguments com- 
pletely independent of the Eddington limit. For example, 
one may demand the value of ^gw at r.m given by (j43|) to 
be equal to some characteristic time t. This would im- 
ply that after the GW emission becomes the dominant 
cause of the binary inspiral, the lifetime of the binary 
until merger is equal to t. Such estimate gives 

7/8 
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Fig. 1. — Properties of a steady disk around a SMBH bi- 
nary with mass Mc = 10'' Mq described by a solution II13I I 
with Fjfi = 5 X 10'*® org and mass accretion rate at infinity 

^J,l = = A'/gjjjj (radiative efficiency of 0.1 is assumed). Solid 
(black) curves describe the run of (a) aspect ratio h/r, (b) tem- 
perature T, (c) surface density E, and (d) angular momentum flux 
Fj with r for a disk with fe = in radiation pressure dominated 
part of the disk. Dashed (blue) curves describe the same but for 
a standard constant M disk with the same mass accretion rate 
Afoo. Dotted extensions of these curves at small r correspond to 
6 = 1 in radiation pressure dominated part of the disk. Open dots 
on each curve correspond to the transition between the radiation 
and gas pressure dominated regimes, while filled dots describe the 
transition between the electron scattering and free-free opacity. 



It sets a lower limit on the value of Fj necessary for 
equal mass {qs = 1) Mc = 10^ Mq SMBH binary to 
merge within the Hubble time. For Fj = i^j.Edd such 
a binary would merge within 4 x 10^ yr after the GW 
emission starts to dominate its orbital evolution. 

Finally, using equations ^ and (pS)) one can also ex- 
press Fj via the disk mass enclosed between its inner 
edge and some outer radius via 



isk 



2Fj 
3a 



Ft 



dr 



(47) 



where it is assumed that gas pressure dominates (or b — 
1), Cs = {kT/ fiy/^) and Fj — const between and r^. 
The approximate relation in this formula is valid for Tq 3> 
Tin , provided that the outer regions of the disk dominate 
its mass — a rather natural assumption as long as the 
midplane temperature falls with increasing radius, as this 
equation shows. Note that Mdisk explicitly depends only 
on the disk temperature at the outer radius but not Tq 
itself. 

Despite the estimates ([^5]) and (|T7)) which have clear 
physical meaning, we still advocate the use of i^j.Edd 
as it represents an important upper limit on Fj: for 
Fj ^ Fj^-Edd the orbital evolution of the binary is essen- 
tially not affected by the Eddington limit (even though 
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at late stages of inspiral the inner part of the disk may 
become super-Eddington and be depleted by the radia- 
tion pressure driven wind, see fJEOiS^USI)- 

3.1.5. Global properties of steady state disks. 

Results of ^ H33.13.1.m33.13.1.3l allow us to under- 
stand global characteristics of steady circumbinary disks 
described by the solution (fT3|) . 

In Figure [1] we demonstrate how the main properties 
of such a disk vary as a function of r across regions with 
different opacity and pressure behavior. We consider a 
Mc = 10^ Mq SMBH binary orbited by a disk described 
by the relation ^ with Fj^o = 5 x 10''^ erg and Fj^i = 

Mao = M^dd- Far from the binary (at r — > cx)) disk 
transitions to a standard constant M = Moo accretion 
disk, while at small separations corresponding to I < 
Fj^o/Moo it becomes a constant Fj disk. Note that M = 
Mao at all radii meaning that mass has to be removed 
at rin at the same rate with which it is supplied at large 
radii. 

We compare the properties of this disk to a standard 
constant M = Moo disk in which Fj^q = (dashed 
curves). We find that the midplane temperature T is 
higher in a disk with non-zero Fj q in all regimes. As a 
result, this disk is more extended vertically and becomes 
geometrically thick (h/r ^ 1) at r = 3 x 10"^ pc (« 300 
Rs) in the radiation pressure dominated regime. Unless 
the disk is truncated on the inside by the binary torque 
at larger separation, it would presumably be losing mass 
at this point. Another consequence of higher T is that 
all transitions between different regimes occur at larger r 
(typically by a factor of several) in the disk with non-zero 
Fj^o- In particular, radiation pressure starts to dominate 
in this disk at « 5 x 10~^ pc as opposed to « 10~^ pc in 
the Fj^o = disk. 

The comparison of the surface density structure be- 
tween the two disks depends on whether 5 = 1 or in 
the radiation pressure dominated regime. In the former 
case S keeps increasing towards small r in both types 
of disks, which reflects the inefficiency of viscosity pro- 
portional to (relatively small) gas pressure in the b — I 
case. We find S to be higher in the Fj^q ^ disk by an 
order of magnitude in the radiation pressure dominated 
regime. 

In the case of v scaling with the radiation pressure 
(6 = 0) the behavior of E is completely different as it falls 
towards small r, because of rapid increase of T and in 
the radiation pressure dominated regime. As a result, 
close to the binary surface density is much smaller (by 
more than an order of magnitude) in the disk with non- 
zero Fj^Q. The relatively small amount of mass residing 
in the vicinity of the SMBH binary in this case may have 
implications for the properties of the afterglow following 
the binary merger. 

It is well known (Goodman 2003) that the outer parts 
of disks around SMBHs may be subject to gravitational 
instability (Safronov 1960; Toomre 1964; Goldreich & 
Lynden-Bell 1965) when the Toomre Q parameter de- 
fined as Q = ^Ics / {ttCS) drops below unity. In the case 
of constant Fj disks or steady disks with Fj{l) given by 
the solution ([T3l) one can easily determine the radius r'^^ 
at which the disk becomes self-gravitating using our re- 
sults for T(r) and S(r) as a function of Fj derived in 



3.1.1113 3.1 3.1.31 for different physical regimes. In 
the interest of brevity we do not perform this straight- 
forward exercise here. For constant M disks described by 
the solution (|13p with Fj q = corresponding expressions 
for r*^ can be found in Goodman (2003) and Haiman et 
al. (2009). Other effects that may invalidate our treat- 
ment at large separations — low optical depth of the disk 
when S becomes small, neutrality of the disk at low T, 
etc. — have been previously discussed in Haiman et al. 
(2009). 

4. DISK EVOLUTION. 

Steady state solutions of the master equation dis- 
cussed so far require rather special circumstances to be 
realized globally, such as M which is independent of r. 
In real circumbinary disks this condition is difficult to 
realize, simply because the mass supply rate of the disk 
at large separations Moo is determined by processes that 
have nothing to do with the central binary. At the same 
time M in the inner disk is set by the binary torque, and 
is in general different from Moo- 

Because of this mismatch of M in different parts of the 
disk gas has to accumulate somewhere and this naturally 
leads to the evolution of disk properties. As we will see 
later evolution typically leads to the establishment of the 
quasi-steady state in the inner parts of the disk, where 
one can then apply the results obtained in previous sec- 
tion. To understand this process we need to obtain the 
time-dependent solutions of equation Q. This neces- 
sarily requires specifying the dependence of the diffusion 
coefficient Dj given by equation ([5|) upon I and Fj and 
a set of boundary conditions. 

With a rare exceptions discussed in !j44.2l we will al- 
ways assume that the circumbinary disk starts out as a 
conventional constant M disk in which 

Fj{t = 0,l)=Mool, (48) 

see equations (fT4|) and ([TTl) . This is a natural assump- 
tion for a viscously evolving disk, which is supplied by a 
source of mass located at r ^ r^. Indeed, as the mass 
flows in, the viscous time in the disk decreases with r 
causing its structure to converge to that of a standard 
constant M disk. When the inner edge of the disk ap- 
proaches the semi-major axis of the binary, the latter 
starts tidally interacting with the disk providing a source 
of angular momentum. Tidal torque on the disk rapidly 
increases and stops the mass inflow at some radius r^ 
(or, equivalently, some value of the speciflc angular mo- 
mentum lin), which is comparable to the semi- major axis 
of the binary. We consider evolution of the disk starting 
from this moment and use an inner boundary condition 
in the form 

dFi 

—f ^M{h„)^xMoo, (49) 

Ol l=lin 

see equation (fT2|) . 

In some cases one can solve equation analytically, 
which is useful for qualitative understanding of the more 
complicated situations. In particular, for Dj — const 
with initial and boundary conditions specifled by equa- 
tions and (Hn]) correspondingly one finds the follow- 
ing solution (with I playing role of spatial coordinate) : 

Fj{t,l)^xMool + il~x)Moo 
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Fig. 2. — Time evolution of the angular momentum flux dis- 
tribution Fj(l,t) (l is the specific angular momentum acting as a 
space-like coordinate) in a disk with constant Dj (see Eq. 1501 ) 
and different boundary conditions imposed at the inner disk radius 
Tin, at which the specific angular momentum is M(li„) = 
(solid lines) and M{lin) = O.bMoo {dotted lines). The initial distri- 
bution of Fj(l,0) = Mccl assumes a standard constant M = Moo 
disk and is shown by a dashed line. Snapshots of Fj at different 
moments of time (labeled on the Figure) are shown. 



I - I, 



X l^in + ^m)erf 



1/2 






j exp 


ADjt 


1 



(50) 



This solution is shown in Figure [2] for different values 
of M{lin) at different moments of time. One can clearly 
see that as the time goes by the influence of the central 
binary extends to larger and larger values of /, implying 
also larger distance from the central binary r. The transi- 
tion from the initial distribution of Fj given by equation 
to Fj strongly affected by the binary torque occurs 
at^ I = Zi,fl ~ {Djt)'/\ 

Note that the solution for M(/i„) = (x = 0) clearly 
exhibits an inner region with Fj — >■ const, in agreement 
with our discussion of M = solutions in Sj2X2l At the 
same time the solution for M{lin) — O.SMqo (x = 0-5) 
develops an inner region with Fj linearly increasing with 
I, again in complete agreement with the solution (jl3l) 
with Fj^i ^ 0. Note that the torque exerted on the disk 
by the binary Fj(lin) is quite different in two cases: it is 
smaller for larger values of x which is natural since higher 
X implies less mass accumulation at the inner edge of the 
disk and less torque exerted by the binary on the disk. 
This means that higher x (higher M{lin)) should result 
in slower inspiral of the binary. 

^ Here we follow the notation of IPP, who called the radius at 
which this transition occurs the "radius of infiuence" rinfl. 



We show next that the main features of disk evolu- 
tion illustrated for the D,j — const remain valid for the 
more general behavior of Dj. Following Filipov (1984), 
Lyubarskij & Shakura (1987), Pringle (1991), IPP, and 
Lipunova & Shakura (2000) we now derive self-similar 
solutions for evolving structure of an externally fed and 
centrally torqued circumbinary disk, allowing for the pos- 
sibility of some mass overflow across the orbit of the sec- 
ondary. Based on these solutions we then outline in §44.21 
a general picture of the non-self-similar disk evolution. 

4.1. Self-similar solutions J or circumbinary disk 
evolution. 

In ^J3] we demonstrated that in a variety of situations 
disk properties can be expressed as power laws of dif- 
ferent physical parameters — Fj, r (or, equivalently, I), 
Mc, etc. Based on that we show in Appendices lAl fc IB] 
that the diffusion coefficient in different regimes can be 
generically expressed in terms of the angular momentum 
flux Fj and specific angular momentum / in the power 
law form 

7d 



Dj = Dj^o F, (51) 

where d and p are constant power law indices. The ex- 
plicit expressions for Dj^, d, and p in relevant regimes 
can be found in Appendix |B] 

Equation (jlj with Dj in the form ([5T|) admits a self- 
similar solution^ provided that the problem at hand has 
no intrinsic scale. In circumbinary disks the inner edge 
of the disk (comparable to the semi-major axis of the 
binary) sets a natural scale. However, after the angular 
momentum injected by the binary has been transmitted 
by viscosity out to distances large compared to the radius 
of the inner cavity, this scale should not affect system's 
behavior and evolution becomes self-similar. 

To illustrate this point let us consider the solution ([50]) 
corresponding to a particular case d ~ p = in the 
hmit of t > Ifn/Dj, so that ^infl > kn (or ri„fl > ri„). 
In this limit we can rewrite the solution (1501) as 



Fj{tJ) = MooVD7t 



xe + (i 



x)Cerf ( I 



2 / ^2 



(52) 



where ^ — l/\/Djt is the dimensionless coordinate, 
which, as we will see later, plays the role of an indepen- 
dent self-similar variable. It is clear from this result that 
at late times the solution of the evolutionary equation ^ 
for Dj = const is indeed independent of the exact value 
of at which the inner boundary condition is imposed, 
and one can effectively set li„ to zero thus eliminating 
any intrinsic scales from the problem at hand. 

In general case of arbitrary d and p in the expression 
(1511) we first define new variables 



fj 



Moo 

transforming equation ^ into 



dr 



fj = P- 



(53) 



(54) 



® Previously Lyubarskij & Shakura (1987) derived such solutions 
for Fj assuming boundary conditions different from what we use 
here. 
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Fig. 3. — Behavior of t he self-similar function / vs the variable 
§ defined in equation I I56I I. Different panels correspond to different 
values of the power law parameters d and p: (a) gas pressure dom- 
inated case with k = Kes, (b) gas pressure dominated case with 
K = K.ff- Various curves represent the run of /(§) for different 
boundary conditions imposed at the inner edge, allowing for the 
possibility of mass inflow past the secondary orbit: M{0) = xMao- 



with the boundary condition 

M{1 



dl 



0) 



i=0 



X, 



(55) 



) , where x ^ 1 ^-^d may be equal 



(instead of equation 
to zero. 

One can easily see that equation ([M)) admits self- 
similar solutions in the form 



1 



d + p — 2 

where function / satisfies ordinary differential equation 

f"fe^n{l-d){f-U') (57) 
with the boundary conditions 

/'(e^o) = x, /'(e^(x)) = i, (58) 

where the second boundary condition follows directly 
from the initial condition (|48p . Note that for d = p = 
equation (j56p reproduces the similarity seen in the ana- 
lytical solution ([52]) obtained for constant Dj. We also 
verified the self-similar scalings of E derived in IPP for 
X = case. 

In Figure ([3]) we show the solutions of this equation for 
different values of mass accretion rate across the orbit of 
the secondary parametrized by the value of x (general- 
izing results of IPP to the case of non-zero x)- We have 
chosen two sets of d and p corresponding to astrophysi- 
cally relevant situations: gas pressure dominated regime 
with K = (d = 2/5, p — —6/5, n — 5/14), and with 
K = Kff {d = 3/10, p — —4/5, n — 2/5), see Appendix IB] 



or Lyubarskij & Shakura (1987). We will see later in fJU 
that the radius of influence often extends into the parts 
of the disk where one of these regimes is valid. 

As expected, the transition from the outer solution 
/(^) ~ ^ unaffected by the binary torque to the inner 
solution influenced by it always occurs at ^ « 1 — 2. This 
transition clearly corresponds to the radius at which the 
viscous time in the disk tytsc is equal to the evolution 
time t of the system, i.e. the time that has passed since 
the central binary started tidally interacting with the 
disk at its inner edge. In other words, tyisdS, ~ 1) i. 

We take the radius of influence rinfl to correspond to 
6nfl = 1- Then, according to equations ([53l) and (|56l) 
the value of the specific angular momentum Una at r-mH 
is given by 



Zinfi(t)= (i^j,oM^)"t", n„fl(i) 



(59) 



while the viscous angular momentum flux in the disk is 

Fj{l, t) = Moohn&it) X / (Z//i„fl(i)) , (60) 

where the value of /(^) can be found from Figure [3] for 
a given x- 

The torque acting on the central binary in the limit 
kn < kn&{t) (or Tin < r'i„fl(i)) is given simply by 
Fj{0,t) — Moohnfi{t)f{0,x) where the dependence of 
/(O, x) upon the efficiency x of mass inflow through the 
secondary orbit is shown in Figure H] for different types of 
disks. One can see in complete analogy with the analyti- 
cal solution ([50)1 that allowing for some mass flow across 
the orbit of the secondary (i.e. assuming non-zero x) 
leads to the reduction of the torque acting on the binary. 

It might seem surprising that the torque experienced 
by the binary is independent of the mass of the secondary 
— Fj{0,t) is set only by Moo and ^infl(i), since it is the 
potential of the secondary that gives rise to the tidal 
coupling with the disk. We comment on this point in Sj5l 

4.2. General description of the disk evolution. 

In H44.ll we outlined main features of the self-similar 
evolution of a circumbinary disk which arises when the 
three essential conditions are met: (1) the behavior of 
the diffusion coefficient Dj is given by a simple power law 
form (|5ip . (2) the outer parts of the disk are well approx- 
imated by a standard constant M disk with Fj = Mool, 
and (3) the radius of influence rinfi far exceeds the semi- 
major axis of the central binary rf, (and the radius of the 
inner disk edge ri„). We now describe how the picture 
of the disk evolution changes when these assumptions 
are relaxed by concentrating on the situation when there 
is no mass inflow across the orbit of the secondary, i.e. 
X = or M{rin) — 0. More complicated setup allowing 
for some mass inflow across the orbit of the secondary 
can be understood by generalization of the picture that 
emerges in Af(ri„) = case. 

Different parts of circumbinary disks can feature dif- 
ferent physical regimes as illustrated in SjS] In this case a 
simple form of Dj given by equation ()5ip will not work 
in the whole disk invalidating global self-similarity. How- 
ever, Dj can still be cast in this form in certain distance 
intervals with power law indices d and p intrinsic to each 
region. For example, one can easily imagine that the in- 
ner disk is in the radiation pressure dominated regime, 
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X=M(0)/M 



Fig. 4. — Value of /(^ = 0, x) ^ function of x — the fraction 
of Moo passing through the orbit of the secondary. Different curves 
correspond to different values of the power lawr parameters d and p 
indicated on panels; {dotted) gas pressure dominated case with k = 
Kes, (solid) gas pressure dominated case with k = kjj. Knowledge 
of /($ = 0) allows computation of the torque acting on the central 
binary for arbitrary value of the "accretion fraction" x- 



while further out it transitions to gas pressure dominated 
regime with the opacity initiahy given by Kes and then 
by K//. In all these regimes one can locally use the power 
law description of Dj as described in Appendix [Bl with 
smooth transitions between the different scalings at the 
boundaries of different regimes. 

Evolution equation written in the form Q easily allows 
us to understand the behavior of the disk properties in 
this more complicated situation. Indeed, interior to rjnfl 
the viscous time in the disk gets shorter with decreasing 
r, meaning that for r < r\-a^ disk tends to approach a 
steady state solution. Then equation ([4]) implies that Fj 
in this part of the disk is given by a simple solution (fT3l) 
independent of the complicated behavior of Dj caused by 
the transitions between different physical regimes. Disk 
properties such as ^(r), r(r), etc. can be computed 
as functions of this radially constant Fj (for M(rm) = 
0) using formulae derived in and will show different 
dependence on r, Fj and other system parameters in 
different regimes. 

At every moment of time the value of the radially con- 
stant (for r < Tinfl) angular momentum flux Fj is obvi- 
ously set by the the disk properties in a particular phys- 
ical regime corresponding to r ^ ?'infl- This regime can 
change in time since Fj steadily increases and both the 
transition radii of different regimes and rinfl vary. Never- 
theless, it is clear that our results for the self-similar disk 
evolution obtained previously should allow one to easily 
understand even more complicated situations, see fJ51 

Second complication arises if the circumbinary disk 
does not start out as a standard constant M disk with 
the initial distribution of the angular momentum flux in 



the form (I48p but is characterized by some more com- 
plicated initial distribution of Fj{l,t = 0). Again, our 
understanding of the self-similar disks allows us to qual- 
itatively characterize disk evolution in this case. The re- 
gion influenced by the binary torque would still expand 
in time with the dependence knfiit) given by an implicit 
relation 

^txDj {Fjihnti, t = 0), /infl) , (61) 

see equation (|4]), in which we explicitly indicated the de- 
pendence of Dj{Fj, I) on Fj and I. Interior to r-mfi angu- 
lar momentum flux is roughly constant with radius and 
equal to Fj{kr,ii{t),t = 0). 

At the same time outside rinfl(i) the disk will still main- 
tain the distribution of Fj close to the initial distribu- 
tion Fj{l,t = 0) since the viscous time there is long 
compared to the system lifetime (which is also equal to 
ifisc(nnfl))- If the initial distribution of Fj exhibits a 
maximum at some radius r^ax then past the moment 
when Tinfl ~ Tmax the circumbinary disk will turn into 
a decretion disk (Pringle 1991) and the mass accumu- 
lated in the central part of the disk will start flowing 
out, driven by continuing injection of the angular mo- 
mentum by the binary (assuming that the binary does 
not merge by that time). 

Finally, initially the radius of influence may not 
strongly exceed the radius of the inner edge of the disk. 
This is true if the evolutionary lifetime of the system has 
not yet exceeded the viscous time at the inner edge of the 
disk. However, at later times the condition rjnfl > ri„ is 
guaranteed to be fulflUed since r-mfi steadily grows while 
both the binary semi-major axis and can only de- 
crease. As a result, at late times the system should in- 
evitably converge to the self-similar mode of evolution 
(see e.g. evolution shown in Figure ^ or its general- 
izations described above for the more complicated situa- 
tions. 

5. IMPLICATIONS FOR SMBH BINARY EVOLUTION. 

We now apply the results obtained in previous sec- 
tions to the coupled evolution of SMBH binaries and 
disks around them. The two processes — orbital evo- 
lution of the binary and evolution of the disk properties 
— must be considered simultaneously because of their 
mutual influence on each other. In exploring this evolu- 
tion one must pay special attention to the nonlocal effect 
of the binary torque on the disk. 

Indeed, according to the results obtained in 21 the 
value of the angular momentum flux Fj carried through 
the disk near the binary (which determines the orbital 
evolution of the binary) is not set locally but is deter- 
mined by the disk properties at the radius Hnfl, which 
is the outermost radius affected by the viscous transport 
of the angular momentum deposited in the disk by the 
binary. This property, often omitted in previous studies 
of the SMBH binary evolution, is very important as we 
show further. 

We emphasize that the torque acting on the binary, 
which according to equation (fTO|) is equal to Fj{rin), is 
independent of the mass of the secondary Mg and is the 
same irrespective of the mass ratio of the binary q (see 
also IFF). Also, the dependence of Fj{rin) on the total 
mass of the binary Mc arises only because Mc determines 
the angular frequency in the disk. At first sight this may 
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Fig. 5. — Characteristic evolution of the torque Fj (labeled on 
the right axis) experienced by the central binary with Mc = Iff 
Mq and (J = 1 as a function of time. Two cases are shown: Fj(t) 
computed self-consistently accounting for the circumbinary disk 
evolution (thick black solid line) and assuming a standard con- 
stant M disk (blue dotted line). In both cases the mass accre- 
tion rate in the disk far from the binary is the same and equal 
to Moo = 10~^ ME(jd- Note the difference between the values of 
Fj computed in two ways for t > 10^ yr. Also shown is the or- 
bital evolution of the binary (thin solid line) computed according 
to the first, self-consistent prescription for Fj{t). Different regimes 
in which the disk can be present are labeled on the plot (red dot- 
ted lines r'^^'^/s^ and r"^/^ show their boundaries) and the gas 
pressure-dominated case with shaded. Dashed curve 

shows the run of r'°fl(t) in different regimes, and is closely related 
to the behavior of self-consistent Fj{t). Eddington limit becomes 
important bel ow the dot -dashed line in the lower right corner of 
the plot, see i|5 5.2 5.2.51 Black square dot marks the transition 
from the disk- to GW-dominated orbital evolution of the binary. 
See text for more details. 



seem strange since the strength of tidal interaction is de- 
termined by the potential of the secondary. However, one 
has to keep in mind that the full torque exerted by the 
binary on the disk is generally found to scale as (Goldre- 
ich & Tremaine 1980; Papaloizou & Lin 1984; Petrovich 
& Rafikov 2012) 

Fj{nn) (X (62) 

irrespective of the precise form of the torque density dis- 
tribution. Here Eg is the disk surface density just outside 
the region where the binary torques are important, and 
A = |ri„ — r^l is the width of the gap — the separation 
between the secondary orbit and the inner edge of the 
disk. Both Fj{rin) and Sq are set globally at the radius 
of influence. 

Equation (|62p shows that a particular value of the bi- 
nary torque Fj(rin) can be obtained not only by chang- 
ing Mg but also by varying the width of the gap A for 
a given Mg, and this is how the disk-binary tidal inter- 
action self-regulates itself to provide a necessary torque 
on the disk. As Fj{rin) varies in time for a fixed Mg the 



width of the gap should also vary. The same is true if one 
varies the mass of the secondary while keeping Fj{rin) 
constant — the width of the gap would simply scale as 
A c>c q^/^ . Of course, q cannot be arbitrarily small since 
a low mass secondary may not be able to prevent the 
mass flow across its orbit (e.g. if A needed to provide a 
given value of Fj{rin) turns out being smaller than the 
disk scaleheight), meaning that the gap does not exist 
in the first place (see ^15 5.2 5.2.6p . But as long as the 
gap opening conditions are satisfied for a given q, the 
width of the gap should always be able to adjust itself to 
provide just the right amount of torque on the disk. 

In our subsequent calculations we will not be directly 
addressing the "final pc" problem (Lodato et al. 2009) 
as we typically follow SMBH binaries starting at rather 
small separations, 10^^ — 10""*. Such binaries may be 
created by previous (possibly multiple) episodes of gas 
infall into the center of the galaxy in which the binary 
resides, each of which would tighten its orbit. At the 
same time some of our findings (e.g. significant reduction 
of the binary inspiral timescale when the disk evolution 
is self-consistently included) are likely to be relevant for 
attempting to resolve the "final pc" problem by account- 
ing for the possibility of a circumbinary disk surrounding 
the binary. 

In all our calculations we take the viscosity to scale 
with the total rather than the gas pressure in the radia- 
tion pressure dominated regime, i.e. 6 = 0. We will also 
assume that the tidal torque of the binary presents suffi- 
ciently strong barrier to infiowing gas to completely sup- 
press gas overflow across the orbit of the secondary. This 
means that the boundary condition at the inner edge of 
the disk is given by x = M{rin)/ M^o = 0. In principle 
one can easily extend our results on the orbital evolution 
of the binary to the case of % 7^ 0. Finally, even though 
the binary itself is not accreting when M(ri„) = the in- 
ner parts of the disk are still gaining mass, which changes 
the potential in which gas orbits further out. In this work 
we are mainly concerned with the disk-related effects on 
the binary evolution and for that reason we neglect the 
increase of the binary-f disk system mass throughout the 
calculation. 

5.1. Binary inspiral: basic features. 

We now look at the details of the orbital evolution 
of the binary SMBH embedded in a circumbinary disk. 
In Figure [S] we show the joint variation of the binary 
and the disk characteristics. The binary orbit is evolved 
according to 



dt taw tj ' 
where ^gw and tj are given by equations (|1T|) and (|^ . 
Variation of the disk properties, including the evolution 
oi Fj{rin) entering the equation (j42|) . is described below. 

We start an equal mass (g = 1) binary with the total 
mass Mc = 10'' Mq with initial semi-major axis of 10"'^ 
pc. At time t = disk properties correspond to a stan- 
dard constant M = MEdd disk extending from very large 
distances (effectively from infinity) down to the binary 
semi-major axis (for simplicity we disregard the differ- 
ence between the binary semi-major axis 7'^ and the inner 
radius of the disk, which is a factor of 2 uncertainly at 
most, see MacFadyen & Milosavljevic 2008). 
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This setup would naturally arise if the binary initially 
resided in a gas-free environment and then gas started 
flowing into the galactic center in a disk-like configu- 
ration from large distances. Because viscous evolution 
accelerates at small radii the disk would naturally set- 
tle into a constant M configuration. At some point its 
inner radius would reach the vicinity of the binary and 
tidal torque would stop the gas inflow. This moment 
represents the starting point for our calculations. 

The radial dependence of disk properties at time t = 
can be found in Haiman et al. (2009) or by setting 
Fj = Moo(GMcr)^/^ as appropriate for a standard con- 
stant M disk in formulae derived in ^'S 3.11 As Figure [5] 
demonstrates the binary starts in the radiation pressure 
dominated part of the disk but the transition to the gas 
pressure dominated regime occurs not too far outside of 
Tf), at r'^'^'^/s^^ ^ 2.5 x 10~^ pc. Opacity switches from 
being dominated by the electron scattering to free-free 
opacity at r'^'^^^ « 0.03 pc. These regimes are clearly 
labeled in Figure [5] 

We also show the run of the radius of influence rinfl(t) 
in time by a dashed curve. At r = rinfi{t) the local vis- 
cous time equals the time since the start of the evolution 
t but this is meaningful only if at time t the disk extends 
to r < rinfi(i), which is not always the case. Neverthe- 
less, this dependence is still a useful concept as it allows 
us to see important transitions in the disk properties if 
we were to take an initial SMBH semi-major axis rb{0) 
different from the value shown in Figure O The actual 
dependence rinfl(t) is calculated using definition ([59]) and 
the expressions for Dj^q, d and n that can be found in 
Appendix [BJ Since these expressions are different in var- 
ious physical regimes the behavior of rinfl(t) exhibits dis- 
tinct transitions as it crosses the boundaries of different 
regimes, clearly visible in Figure [3 

As the binary starts tidally interacting with the disk 
the infiowing material begins accumulating near the in- 
ner edge of the disk. In addition to the evolution of 
the binary orbit (thin solid curve) this Figure also shows 
the time behavior of the torque Fj{rin,t) at the in- 
ner edge of the disk (thick solid curve calibrated on 
the right axis), which is absorbed by the binary and 
causes its orbital evolution. The Fj{rin, t) curve initially 
closely follows'' rb{t) because the torque exerted on the 
binary by a constant M disk is well approximated by 
Fjijimt) = MoaiGMcTbY^'^ as long as the inner edge 
of the disk tracks the binary orbit. This torque is small 
enough for rf, not to change appreciably for rather long 
time. 

Approximately at = 2 x 10*^ yr, when i^isc(''&(0)) ^ t, 
the radius of infiuence of the binary torque rjnfl grows 
beyond the initial binary semi- major axis rf,(0). As a 
result, angular momentum flux at the inner edge of the 
disk Fj(rin,t) starts increasing as Fj{rin,t) cc imdit) cc 
t^/"^ thus accelerating the orbital evolution of the binary. 
Unlike the calculations of the self-similar disk behavior 
in ^ 14 4.1l in our present calculation we do not consider the 
details of the smooth transition between the parts of the 
disk inside and outside of rinfl. Instead we simply assume 
that Fj{r,t) — Mool for I > /infl(i) as appropriate for a 

For clarity we have slightly shifted the curve Fj {rin , t) upward 
in Figure[5]to avoid overlap with other curves, e.g. ri,{t). 



constant AI disk, while Fj{r,t) = Moohnfi{t) is constant 
in space for I < Zinfl(i). In other words, we adopt a simple 
piece- wise dependence of Fj on r: 

FAr,t)^l M^iGM^rf\^ r > n^t), (g^) 

For that reason at late times Fjijimt) starts tracking 
the run of rinfi(t) in Figure [H 

Initially rjnfl stays in the radiation pressure dominated 
regime and varies as rinfl oc t^/"^ , see equations (f59| and 
(jBip . At t « 2 X 10^ yr the radius of influence reaches 
out into the gas pressure dominated regime with k = Kes ■ 
There Hnfl grows as rinfl oc t^^"^ (see the break in slope 
of Tinfl curve), and the increase of Fj{rin,t) and of the 
torque acting on the binary accelerates. 

In our calculations the disk interior to r = rjnfl is a 
constant Fj disk (see equation [64]) with properties ex- 
plored in 3.11 rather than a constant M disk present 
outside of the radius of influence. As the value of Fj 
in the inner disk grows the disk properties keep chang- 
ing as well. In particular the boundaries of the different 
regimes, i.e. ^''^^'^/sas g^^^j j.es/ff ^ expand as the time goes 
by. Thus, the state of the disk at the location of the 
binary may change not only because of the variation of 
the binary orbit, but also due to the disk evolution, see 
Figures [6]& in 

Note, that even though in the calculation shown in Fig- 
ure [S] the binary is always in contact with the radiation 
pressure dominated region of the disk, beyond i « 2 x 10^ 
yr the torque on the binary is determined by the gas pres- 
sure dominated part of the disk with k = Kes , since this 
is where rjufl is. This demonstrates the nonlocality of 
the disk influence on the binary — the state of the disk 
near the binary is essentially irrelevant for its orbital evo- 
lution. It is what goes on in the disk at r ~ rjnfl that 
determines the torque acting on the SMBH binary at late 
times. 

This fact has been overlooked in previous studies of 
the gas-assisted SMBH inspiral problem. In particular, 
Haiman et al. (2009) and Kocsis et al. (2011) used the 
self-similar results of SC95 and IFF to account for the 
mass accumulation in the disk near the binary orbit. 
However, in their calculations they have effectively as- 
sumed that the radius of influence rinfl corresponds to 
the same regime of the disk in which the binary is cur- 
rently residing. In the case shown in Figure [5] this would 
mean that even after 2 x 10^ yr the increase of rjnfl with 
time would be calculated according to the scaling for the 
radiation pressure dominated regime, leading to an un- 
derestimate of the torque acting on the binary and an 
overestimate of its inspiral time. Our calculations fully 
take into account the nonlocality of the disk influence on 
the disk. 

To better illustrate the role of mass puleup and disk 
evolution for the orbital evolution of the binary we com- 
pare our results with calculations in which disk proper- 
ties remain well represented by the properties of a con- 
stant M disk at all times. Then the torque on the bi- 

1/2 

nary is always given by Fj{t) = M^o [GMcri,{t)] , see 
dotted curve in Figure [3 Clearly, this torque can be 
much smaller than the real Fj{t), especially at late times, 
meaning that such calculations should overestimate the 
inspiral time of the binary. 
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Fig. 6. — Comparison of the orbital evolution of the SMBH binary computed in two ways: fully accounting for the disk evolution driven 
by the tidal torque of the binary (solid black curves) and simply assuming the disk properties to be given by a standard constant M solution 
(dotted blue curves). In both cases the value of Moo, initial semi-major axis of the binary, and all other parameters are assumed to be the 

same. Results are shown for Mc = 10^ Mq, Mac = A^Edd (panel a) and for Mc = 10^ Mq, Moo = M^^^ (panel b), with binary mass 
ratio q = 1 in all cases. Meaning of other curves and labels on these plots is the same as in Figure [5] Note that properly accounting for 
the binary-driven disk evolution can reduce the lifetime of the binary by almost an order of magnitude in some cases. 



In Figure [6] we display ri,{t) calculated using constant 
M disk properties (dotted curves) and fully accounting 
for the binary-driven disk evolution (solid curves) for two 
equal mass {q — 1) SMBH binary+disk systems: one 
with Mc = 10^ M0, Moo = MEdd and another with 
Mc = 10"^ M0, Moo = 10"^AfEdd- Two different starting 
radii rf,(0) are explored in both cases. 

One clearly sees that in most cases evolution of ri, com- 
puted using standard constant M disk properties consid- 
erably overestimates (by almost an order of magnitude 
in some cases) the binary inspiral time. In the case of 

Mc = 10^ Mq, Moo = lO-^A/Edd, and n{0) = lO'^ pc 
the evolutionary tracks of the binary computed by two 
methods coincide. This is because in this case GW emis- 
sion is more important for the evolution of the binary 
than the tidal coupling to the disk at all times, result- 
ing in a universal behavior of rb{t). But whenever disk 
torques are important, binaries shrink faster when the 
self-consistent disk evolution is properly taken into ac- 
count. For that reason, we strongly discourage the use 
of standard constant M disk properties for exploring the 
evolution of the central binary of a circumbinary disk. 

5.2. Binary inspiral: parameter exploration. 

In Figures [7] & [8] we provide a more detailed and sys- 
tematic view of the SMBH binary evolution under differ- 
ent conditions. We explore two representative values of 
Mc- 10^ Mq (implying Schwarzschild radius Rs = 10~* 
pc) and 10^ Mq [Rs = lO'^ pc), but our results can 
be trivially extended to other values of Mc. Accretion 
rate in the disk at large distances is taken to be either 



Moo = lO^^A/Edd or MEdd, and binary mass ratio is 
varied between q = 10~^ and 1. We also consider two 
different values of the starting semi-major axis of the bi- 
nary rbiO): 10-* pc and 10"^ pc for Mc = 10^ Mq, and 
10-3 pc and 10"^ pc for Mc = 10^ Mq. These values 
are close to the "bottleneck" semi-major axes at which 
the stellar dynamical orbital evolution of SMBH pairs 
decelerates dramatically, see Yu (2002). 

One might worry that the outer parts of the disk can 
be prone to gravitational instability. We determined that 
in the initial constant M disk Toomre Q equals unity at 
r'^f = 9 X 10"^ pc, 0.3 pc, 0.3 pc, and 0.5 pc for the sys- 
tems shown in panels (a)-(d) correspondingly in Figure 
[7] (for a = 0.1, e = 0.1, = 0.5mp). Thus, at the start 
of our calculations one needs to worry about the impor- 
tance of the disk self-gravity only for the system with 
Mc — lO"^ Mq, Moo = M^dd and the binary starting at 
rb{0) = 10~^ pc, where the disk can be marginally grav- 
itationally unstable. One has to keep in mind though 
that later the torque on the binary is going to be set at 
the radius of influence that expands beyond ^(,(0), see 
^5 5.11 If Tinfl would exceed r^^ at some point, the calcu- 
lations of disk evolution would need to be refined. In our 
present study we neglect this complication; it may only 
be an issue for the high-mass systems shown in Figure 

Ell. 

Figure [7] shows evolutionary tracks of the binary orbit 
mapped onto the disk state in the format analogous to 
that used in Figures [S] & HI In addition, in Figure [B^,b 
we show the dependence of the binary orbital evolution 
timescale tev = —rb/ri, vs. the binary orbital period Porb 
or semi-major axis rf, (upper axis) for the evolutionary 
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Fig. 7. — Evolutionary tracks of SMBH binaries for different values of Mc = , 10^ M© and Moo/MEdd = 1, 10~^, labeled on 
panels. Tracks for g = 1 (solid black) and 10~^ {dotted red) are shown and in each panel we consider two starting values of the binary 
semi-major axis r^,{0). Layout of these plots and meanin g o f different curves and labels are the same as in Figure|5] In particular, dashed 
and dot-dashed curves show the run of ring (equation [59) ) and rgdd (equation [65]) correspondingly, and blac k square do ts mark the 
transition to the GW-dominated orbital decay. Open pentagons mark the possible onset of the gas overflow, see ^5 5.2 5.2.61 See text for 
detailed description of evolution. 



tracks displayed in Figure [71 This plot allows us to eas- 
ily see the transition from the disk dominated evolution 
at longer periods to the GW dominated phase, which is 
clearly described by straight line tracks at small values 
of Porb- Figure [8K,b can be directly compared to the 
analogous tev{Porb) plots in Haiman et al. (2009). 

Results presented in Figures [7] & [8] can be summarized 
in the following set of conclusions. 

5.2.1. Circumbinary disk can be efficient in driving 
orbital evolution of the binary. 

We generally agree with the results of existing studies 
(e.g. IPP, Haiman et al. 2009) that disks can appreciably 
accelerate orbital evolution of SMBH binaries. Figures 



[7| & [HI clearly illustrate this point. According to equa- 
tion (|4T|) without the disk the orbital evolution timescale 
due to emission of gravitational waves is rather long for 
some of the systems shown in these plots: for example, 
^GW = 8.3 X IQio yr for equal mass 10^ Mq SMBH bi- 
nary starting at r&(0) = 10~^ pc. Systems with small 
mass ratios evolve even slower: a binary with q — 10~^ 
takes ~ 25 times longer to merge due to the GW emission 
alone. The only system in our sample strongly affected 
by the GW emission from the very start is the equal mass 
10'' M0 binary starting at 10"'^ pc and surrounded by a 

disk with Moo = 10~^AfEdd, see Figure [8^. As expected 
this system merges faster than its q ~ 10~^ counter- 
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Fig. 8. — (a,b) Evolution timescale tev = \dlnri,/dt\~^ of SMBH binaries as a function of orbital period Parb or semi-major axis r;, 
(upper axis; use to read off the initial semi-major axis ri,(0)). Starting points for each evolutionary track are marked by either triangles 
(for AfcxD = 10~'^MEdd) Of hexagons (for Moo = ^Edd)- Solid (black) and dotted (red) curves are used for 9 = 1 and q = 10~^ binary 
tracks. Straight portions of these tracks at low P^rb correspond to the epoch when the orbital evolution of the binary is dominated by 
the GW emission. Note that prior to that tracks starting at different rf,(0) do not overlap, illustrating hysteresis of the binary evolution 
caused by the non-local nature of the self-consistent disk-binary coupling. (c,d) Ratio of the local disk mass to the secondary mass Ms 
for a sub-sample of evolutionary tracks displayed in Figure[7l showing that M^/Ms 1 in our calculations. Different tracks are labeled as 

( logtA^oo/A/Eddli logq, log[r(,(0)/pc] J . See text for details. 



part, even though the latter is affected more by the disk 
torques. 

On the other hand, the same q — 1, 10^ M© binary 
starting at rb(0) = 10~^ pc and surrounded a disk ac- 
creting at Moo = -MEdd merges within 7 Myr — more 
than 4 orders of magnitude faster than without the disk! 



Lower mass disks are of course less efficient at driving or- 
bital evolution of SMBH binaries — the same binary sur- 
rounded by a disk with M^o — 10~^MEdd merges within 
3 X 10* yr, but this is still much shorter than the cor- 
responding ^GW- Lower mass ratio binaries are affected 
by the disk even stronger, provided that they present an 
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efficient barrier to the mass inflow at the orbit of the 
secondary. For example, lO*" Mq binary with q = 10^^ 
surrounded by a disk with Moo — 10~^AfEdd, merges 
within 3 x 10^ yr, about an order of magnitude faster 
than the q — I binary with the same parameters. 

In Figure [7] black square dots mark the location on the 
evolutionary track of each SMBH binary where the disk 
dominated evolution switches to the GW dominated or- 
bital decay. In systems with massive, high- Moo disks this 
transition typically occurs when the binary is in contact 
with the radiation pressure dominated part of the disk. 
However, in systems with less massive, lower- Moo disks 
this transition may happen even while the binary is sur- 
rounded by the gas pressure dominated (with k = Kes) 
part of the disk. 

Another way to state the importance of the disk dom- 
inated evolution is to note that in most cases transition 
to GW dominated regime occurs at ri, ^ ''^(O), i.e. af- 
ter the binary semi-major axis has been significantly re- 
duced by the disk torques. For example, evolution of a 
Mc = 10^ Mq binary with q = 10^^ starting at 10~^ 
pc is dominated by torques produced by a A/qo = -^^Edd 
disk down to r?, 2 x 10~^ pc, see Figure [T]:. This is 
almost three orders of magnitude smaller than rt, (0) and 
is about 200 Rs for the binary. This additionally empha- 
sizes the important role of disk torques in shrinking the 
SMBH binary orbits, even at relatively large separations. 

5.2.2. Nonlocal character of the disk-binary coupling 

As described in §5 5. li the torque exerted on the binary 
is set by the disk properties at rinfl. Initially binary can 
affect only its immediate surroundings as it takes certain 
time for the disk to absorb enough angular momentum 
injected by the binary to affect the surface density distri- 
bution further out in the disk. For that reason initially 
Fj{rin) is essentially the same as in the case of a con- 
stant M = Moo disk, and is set by the disk locally, at 
r ^ r^. This torque is usually rather small implying slow 
orbital evolution and long values of t^v ■ 

However, after system has evolved for time comparable 
to the viscous timescale at the inner disk edge tyisdrin), 
Tinfl starts exceeding rf,. Past that point the torque on 
the binary Fj{rin) is being set globally, at distances far 
exceeding r^, which can be clearly seen in several evo- 
lutionary tracks shown in Figure [T] Increase of Fj{rin) 
initially occurs at almost constant ri,, which is reflected 
in almost vertical initial evolutionary tracks in t^v ~ Porb 
plane shown in Figure [8K,b, clearly visible for Mc — 10^ 
Mq and low Moo/MEdd = 10"^. 

Global nature of the torques is generally more pro- 
nounced for lower mass binaries and for lower Moo. This 
is because higher Mc implies earlier transition to the 
GW-dominated orbital decay (see Figure |8K,b) , short- 
ening the binary lifetime and not allowing rinfl to extend 
as far as in the lower Mc case. Higher MooZ-^fedd plays 
similar role, shortening the binary lifetime and reducing 
the radius of influence at the end of inspiral compared to 
the lower Moo case. This is true even though rinfl itself 
depends on both Mc and Afoo — these dependencies are 
usually rather weak, see Appendix |B] 

Almost all tracks shown in Figure [7] at some point in 
their evolution run into the situation described in !j5 5.1[ 



where the inner edge of the disk (assumed equal to r?,) 
and Tinfl reside in parts of the disk corresponding to dif- 
ferent physical regimes. This means that calculations of 
the binary evolution assuming the scaling of rinfl(0 
always correspond to the physical state of the disk near 
the binary orbit (Haiman et al. 2009; Kocsis et al. 2011) 
are not accurate. Figure [7] demonstrates that using this 
simple-minded procedure for a binary in the radiation 
pressure-dominated regime can easily underestimate the 
torque Fj{ri„) driving its orbital evolution thus overes- 
timating the binary lifetime (see e.g. Figure [7^, c). Anal- 
ogously, for a binary residing in the gas pressure domi- 
nated regime with k — Kes and rinfl extending into the 
part of the disk where k = Kff (see Figure [7}D-d) the use 
of nnfi(^) scaling corresponding to the binary location 
would again lead to a (mild) underestimate of Fj{rin) 
(nnfl(0 ^^^^ when K = Kff) and an overestimate of the 
merger time. 

Another consequence of the non-locality of the disk- 
binary coupling is the clear hysteresis in the evolution 
of the system — the dependence of the current rate of 
orbital decay of the binary on the previous history of 
the disk evolution. This property is most readily seen 
in Figure [H in which the evolutionary tracks computed 
for the same Mc, q, M^o but starting at different initial 
radii do not result in the same orbital decay timescale tcv 
at a given orbital period in the disk-dominated regime. 
One can see that for a fixed Porb the inspiral timescales 
computed for different starting conditions can differ by a 
factor of several. This is in contrast to local calculations 
presented in Haiman et al. (2009), in which orbital evo- 
lution depends only on the current value of rf,, see their 
Figures 1-5 showing just a single evolutionary track for 
a given set of Mc and q. 

5.2.3. Evolution of the disk accelerates orbital 
evolution of the binary 

Because the radius rinfi setting the value of the in- 
ner torque Fj{rin) in our calculations (see equation [64] ') 
steadily grows, we necessarily find the disk evolution to 
result in the speed up of the binary decay compared to 
the case in which the disk is a constant M = Moo disk 
at all times. As we have shown in §5 5.11 this results in 
shorter lifetime of the system in our calculations, which 
is a very natural result. 

Interestingly, previously Haiman et al. (2009) have 
reached a directly opposite conclusion — that the disk 
evolution caused by the mass pile up near the orbit of 
the secondary slows down its inspiral (see their Figures 
6 & 7). Whether this difference is caused by the local 
character of the of the disk-binary coupling assumed by 
Haiman et al. (2009) or by the adoption of SC95 solution 
in their calculations is not clear. 

5.2.4. Validity of the secondary dominated regime 

In Figure [8}:,d we display the ratio of the secondary 
mass AIs to the local disk mass Md = S(rh)r^. Our 
calculations explicitly assume a "secondary-dominated" 
hmit Md/Ms < 1 (SC95; IFF; Haiman et al. 2009) since 
only in this case the boundary condition in the form 
(fTO|) is valid, as discussed in §2 2.11 One can easily see 
that almost all evolutionary tracks shown in Figure [5J;,d 
(which are also present in Figure [7]) satisfy Md/Ms ^ 1 
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and thus correspond to the secondary-dominated regime. 
This finding is in agreement with Haiman et al. (2009) 
conclusion on the ubiquity and importance of this stage 
of the binary evolution. 

Evolutionary track for a g = 10~^, Mc = 10 Mq 
SMBH binary starting at 10~^ pc and surrounded by 
a disk accreting at Moo = -A^Edd shows the highest 
Md/Ms ~ 0.2 at early stages of the evolution among 
all the tracks computed in this work. For this binary the 
secondary-dominated regime is close to being marginally 
violated, and its track exhibits noticeable evolution even 
prior to crossing rinfl(i) curve. Whenever ~ Mg the 
inner part of the disk contains the amount of angular mo- 
mentum comparable to the full orbital angular momen- 
tum of the binary at the very start of evolution. Thus, 
to cause appreciable evolution of the disk via tidal cou- 
pling (which happens on the local viscous timescale) and 
to push the radius of influence beyond rh(0) requires the 
binary to give off a noticeable fraction of its orbital an- 
gular momentum. As a result a massive disk is able to 
shrink the binary orbit very efficiently, in less than the 
local viscous timescale at the inner edge of the disk. 

At even higher values of Md/Mg one needs to employ 
a boundary condition ([7]) different from (|10p to properly 
compute the disk evolution. 

5.2.5. Importance of the Eddington limit 

As we discussed in ^ 33 3.1 3.1.41 the magnitude of the 
angular momentum flux in a constant Fj disk cannot 
be arbitrarily high, because of the existence of the Ed- 
dington limit ()40p in the radiation pressure dominated 
region of the disk. We now check how important is this 
constraint for the SMBH binary evolution. 

In Figure [7] we plot as a dot-dashed curve the Edding- 
ton radius rEdd at which the Eddington limit becomes 
important. This dependence can be easily derived from 
equation (j40|l and reads 
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where the value of Fj in the inner part of the disk af- 
fected by the binary torque is assumed to be given by 
equation (|64| . In this expression we use the dependence 
nnfl(i), which is also plotted in the same Figure. Edding- 
ton limit is important for a circumbinary disk whenever 

Figure [7] demonstrates that for our choices of Mc, q, 
Moo and rb{0) the Eddington limit is essentially irrele- 
vant during the disk-dominated phase of the orbital evo- 
lution of the binary: dot-dashed curve of rEdd{t) always 
passes below or close to the square dots marking the 
transition to the GW-dominated regime. This provides 
justification for our calculations of the binary evolution 
since the details of the disk physics, including the advent 
of the Eddington limit, are not going to affect it in the 
GW-dominated regime. 

On the other hand, for all evolutionary tracks depicted 
in Figure [7] the Eddington limit does become important 
at some point and affects the disk properties right before 



the binary merger. This may have important implica- 
tions for the electromagnetic precursor of the merger. 

It is worth stressing that the Eddington limit in a con- 
stant Fj disk can be easily reached even if the mass ac- 
cretion rate in the disk at large separations is far less 
than MEdd, see Figure [7}3,d. This is obviously caused by 
accumulation of mass at the inner edge of the disk result- 
ing in higher temperature and more important role of the 
radiation pressure than in a constant M disk. Thus, the 
Eddington limit in self-consistently evolved circumbinary 
disks is more stringent than in the standard constant M 
disks. 

At the same time, the Eddington limit is still less 
important for low values of Moo/ M^dd- for Moo = 
lO^^A/Edd it kicks in only at the distances of several 
tens of Schwarzschild radii at most, see [7lD,d. At the 
same time for Moo = M^dd the Eddington limit starts 
affecting the disk near the binary as soon as the latter 
enters the GW-dominated regime. 

5.2.6. Gas overflow across the orbit of the secondary 

Our assumption* of Mijin) = adopted throughout 
most of this work is equivalent to demanding the width 
of the gap between the orbit of the secondary and the 
inner edge of the disk A to be larger than the disk scale- 
height h. Indeed, the torque density produced by the 
planet drops for |r — rf,| < h due to the phenomenon of 
the "torque cutoff" (Goldreich & Tremaine 1980), which 
implies that the secondary can effectively repel the disk 
fluid only if the gap width satisfles A > h. When this 
condition is not fulfilled, gas enters the torque cutoff zone 
near the orbit of the secondary where the tidal repulsion 
is no longer effective and starts to overfiow the orbit of 
the secondary (Kocsis et al. 2012). As a result, the inner 
boundary condition in the form M(ri„) <C Moo may get 
violated. 

Following SC95 we can estimate the gap width as 
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(66) 



where we assumed viscosity to be given by equation (|23p . 
In both the gas pressure and the radiation pressure dom- 
inated regime with 6 = (which we adopt in our calcula- 
tions) one obtains the same value of A as in SC95. How- 
ever, in the radiation pressure dominated regime with 
h — 1 (as adopted e.g. by Kocsis et al. 2012a, b) the gap 
is wider by a factor of /3~i/'^, making overfiow less likely. 

Liu & Shapiro (2010) calculated mass accretion rate 
across the orbit of the secondary using a local steady 
state model for the disk structure. They found that 
M{rin) is exponentially small when the factor^ in square 
brackets in equation is large (> 5 — 10). On the 
contrary, when this factor is < 5 the mass accretion rate 
across the orbit of the secondary is found to be close to 
Moo- This provides justification for our use of the bound- 
ary condition M{rin) = (x = 0) whenever A//i > 1, 
and for assuming overflow to occur for A//i < 1. 

* I am grateful to Zoltan Haiman for suggesting the overflow 
calculations presented in this section. 

^ Liu & Shapiro (2010) call this factor g. 
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The condition (|66p allows us to find the value of r at 
which A/h — 1, which we call the "overflow" radius rof- 
Assuming that when the overflow begins the disk around 
the binary is in the radiation pressure dominated regime 
(this is always the case in our calculations, see Figure [71) 
with 6 = and h/r is given by equation ()29|) . we find 
that A//i = 1 at 

/ ^ ^ 2/15 

rof ~ r-Edd i-^j , (67) 

where r-Edd is given by equation (p5)) . 

It is clear from this expression that for equal mass bi- 
naries overflow occurs only after the Eddington limit be- 
comes important since then Tq/ < fEdd- On the other 
hand, as long as g < a^/^ one finds rof ^ fEdd but the 
actual value of rof never deviates too much (i.e. not by 
orders of magnitude) from rsdd because of the weak de- 
pendence of rof on q and a. Indeed, for a = 0.1 and 
q = 10~^ one finds rof ~ 2.5rEdd- 

In Figure [7] we show the position of Tq/ for each evo- 
lutionary track with open pentagons. These tags always 
lie close to the dot-dashed curve showing the Edding- 
ton limit (|65|) . Another important observation that can 
be made by inspecting this Figure is that at least for 
q > 10~^ overflow always occurs close to the start of 
the GW-dominated phase of the orbital evolution of the 
binary. Indeed, only for q = 10^^ and only in massive 
disks with Moo = -^^Edd do we flnd that overflow pre- 
cedes (by only a factor of ~ 2 in terms of rb) the stage 
of the GW-driven evolution, see Figure [7^,c. Both for 
g ~ 1 and for Moo ^ -^Edd (essentially irrespective of q) 
overflow occurs when the orbital evolution of the binary 
is already fully determined by the GW emission. 

Kocsis et al. (2012) have derived quasi-steady solutions 
for the disk structure in presence of the overflow, reminis- 
cent of the SC95 results. For these solutions to become 
valid after the overflow begins, information on the new 
boundary condition at the inner edge of the disk must 
propagate to the current radius of influence r-mfi where 
the unperturbed, standard constant M disk starts. Oth- 
erwise the solution would not converge to a standard con- 
stant M = Moo solution at large radii and M at the inner 
edge of the disk cannot be assumed equal to Moo- Es- 
tablishing connection to the outer disk takes of order the 
viscous time at rjnfl, which is about the age of the system 
in our calculations. Using Figures [7] and [5] one can easily 
see that unless q < 10~^ and Moo — ^^Edd the orbit of 
the binary evolves on a much shorter timescale than the 
viscous time at r-mfi- 

For example, a disk with Moo = ^Edd around a bi- 
nary with Mc = 10"^ M©, q = 10"^ starting at 0.01 pc 
begins to overflow when its period is « 0.1 yr (see Figure 
[7^) and the orbital evolution time scale is « 7 x 10^ yr 
(Figure [Htt) ■ This is almost the same as the viscous time 
at Tinfl for the corresponding evolutionary track. Thus, a 
quasi-steady solution can be marginally valid in this case. 
But if we now look at q = 1 binary keeping everything 
else the same we find overflow to occur at Porb = 0.03 
yr. At this period tev ~ 10'^ yr which is much shorter 
than the viscous time at nnfl (~ 10^ yr). As a result, a 
global quasi-steady solution does not get established in 
this case. Similar situation occurs for Moo ^ MEdd (and 
arbitrary q)- 



To summarize, the overflow across the orbit of the sec- 
ondary is most important for low q, high Moo systems. 
Even then it does not strongly affect the orbital evolu- 
tion of the binary during the disk-driven stage and be- 
comes truly important only when the binary inspiral is 
dominated by the GW emission. These conclusions are 
reached in a setup most favorable for the emergence of 
the overflow — radiation pressure dominated disk with 
v proportional to the radiation pressure pr- If instead v 
scales with gas pressure the overflow is going to be even 
less important for the gas-assisted SMBH binary evolu- 
tion. 

5.2.7. Implications for the gravitational wave 
signatures of SMBH binary merger. 

Our results have interesting implications for future 
space-based gravitational wave antennae such as LISA. If 
we adopt 0.03 mHz (sa lO"'^ yr"^) as a characteristic low- 
est frequency probed by such experiments then according 
to Figure |H^,b equal mass binaries with Mc > 10'' Mq 
are detectable only when their orbital decay is already 
fully dominated by the GW emission, even for massive 

disks with Moo — MEdd- 

However, lower mass ratio systems including the so- 
called extreme mass ratio inspirals (EMRIs) can enter the 
detection band of the GW experiments during the stage 
when their orbital evolution is still dominated by the 
tidal coupling to the disk. An example of this can be seen 
in Figure [Hb where the SMBH binaries with q = 10~^ 
are pushed by the disk with Moo = A^Edd all the way 
until the transition to the GW-dominated regime occurs 
at Porb ~ 10~^ yr. Decay of systems with even lower 
q or Mc can be dominated by their disks down to even 
shorter orbital periods, making detection of disk-driven 
migration quite plausible for low Mc EMRIs (provided 
that the strain they produce is above the signal-to-noise 
of the G W antenna) . 

Disk effects on the GW signal manifest themselves via 
the orbital phase shift of the binary caused by the varia- 
tion of its semi-major axis due to the disk-driven migra- 
tion, see Kocsis et al. (2011), Yunes et al. (2011). Even if 
signatures of the disk-driven migration are indeed found 
in the GW signal of coalescing binaries, it is unlikely that 
one would be able to use these measurements to probe 
the properties of the radiation pressure dominated part 
of the disk in the immediate vicinity of the binary. The 
reason for that again lies in the nonlocal nature of the 
torque acting on the binary: close to Porb = yr 
evolution of all low- Mc binaries shown in Figure [71;,d is 
typically dominated by torques set at rinfl located in the 
gas pressure dominated regime. Thus, GW phase signal 
will contain information only about the properties of the 
gas pressure dominated part of the disk, and will not in- 
form us on the physics of the inner, radiation pressure 
dominated regions. 

We also note in this regard that the calculations of GW 
shifts presented in Kocsis et al. (2011) & Yunes et al. 
(2011) should be revised to account for this non- locality 
of the disk torques on the binary. 

5.3. Spectra of disks around SMBH binaries- 
Disks around SMBHs exhibit a set of observational sig- 
natures which distinguish them from the regular constant 
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Fig. 9. — Evolution of the spectrum of a circumbinary disk 
evolving together with the central SMBH binary of mass Mc = 10^ 
M0 and 5 = 1, starting with initial semi-major axis 10~^ pc. Disk 
has Moo = 10~^ MEdd fa-r from the binary and extends out to 
0.1 pc. Its inner radius evolves together with the binary whose 
evolutionary track is shown in Figure [7Ji. Different curves show 
disk spectrum at different times labeled on the panel. Thick lines 
correspond to a disk self-consistently evolving under the action 
of the binary torque with no mass flow across the orbit of the 
secondary allowed (x = 0). Thin (blue) lines correspond to a 
constant M = Moo disk occupying the same region of space. See 
text for details. 



M accretion disks. We illustrate this difference in Figure 
[S]by showing the spectral energy distribution (SED) of a 
circumbinary disk at different stages of the disk+binary 
evolution. This particular calculation assumes an equal 
mass, Mc = 10^ Mq binary starting at n{0) = 10^* 
pc, surrounded by a disk accreting gas at large separa- 
tions at the rate Moo — lO^^AfEdd; see Figure [TJi for an 
evolutionary track of this system. 

One type of spectra shown by thick curves in Figure 
[3] assumes radial distribution of the angular momentum 
flux in the disk Fj{r, t) to be given by equation (IMl) with 
nnfl(i) taken from the self-consistent calculations pre- 
sented in Figure [Tli. Another set of SEDs (thin curves) 
is computed for the same moments of time (and for the 
same values of ri,{t)) for a standard constant M = Moo 
disk extending down to the binary orbit rf,(<). In this case 

Fj(r,t) = Moo (GMcr)^^^ through the entire disk. Both 
kinds of calculations assume an outer edge of the disk to 
lie at 0.1 pc (in this calculation we disregard complica- 
tions arising at large separations, which are mentioned 
in ^3 3.1 3.1.51 the choice of the outer radius is not im- 
portant to us). We now go over the details of these cal- 
culations. 

After t = 10^ yr of evolution the binary semi-major 
axis has essentially not changed and the radius of its tidal 
influence has extended only out to 3 x 10^** pc. Since this 
value of Tinfl is close to rt, binary torques affect only the 
very innermost part of the disk and the spectra computed 



in two ways (thick and thin long-dashed curves in Figure 
^ do not show significant difference. At A = 10 — 500 
lim both are well fit by a power law vF^, cx v^^^ typical 
for a constant M disk. 

At t = 10^ yr rt is still very close to rb(0), but the ef- 
fects of the binary torque have been viscously transmit- 
ted through the disk out to Hnfl ~ 1.7 x 10^"^ pc. This 
results in a factor of « 4 difference in the torque act- 
ing on the binary in two cases, and noticeably changes 
the spectrum of the disk in a self-consistent calculation: 
the peak wavelength of the spectrum shifts to a slightly 
shorter wavelength and the peak amplitude of vF^, in- 
creases by a factor of 2 compared to t = 10^ yr. Note 
that the latter is a consequence only of the change in the 
disk structure due to binary torques — the inner disk 
radius stays essentially the same in these two epochs. 
For the same reason there is not difference in the spectra 
computed assuming constant M disk (the thin curves for 
t = 10^ and 10^ essentially overlap in Figure E]). 

At t = 10'' yr rt has shrunk to 8.7 x 10~^ pc, while 
Hnfl ~ 10~^ pc. Prior to this moment of time the disk in 
the immediate vicinity of the binary was in a gas pressure 
dominated state with free-free opacity, and the radius of 
influence was also in the same regime. At t = 10^ yr 
binary starts entering the gas pressure dominated part of 
the disk with Figure[7jl. While the spectrum 

of a constant M disk is almost the same at this epoch, 
the SED of a self-consistently evolved disk (thick dotted 
curve) exhibits not only an increase in amplitude and a 
shift towards shorter wavelengths, but also a change in 
slope at A = 10 — 500 ^m: in this range i/F^ is clearly 
steeper than i^'*/'^. All that is again predominantly due 
to the evolution of the radial structure of the disk under 
the action of the binary torque. 

Finally, at t w 4 x 10^ yr binary orbit shrinks to 10^^ 
pc, while Tinfl ~ 0.03 pc. Note that at this separation or- 
bital evolution of the binary is dominated by the GW 
emission rather than the disk torques, see Figure [Tti. 
However, the binary still causes the disk to evolve as 
long as its torque prevents the mass inflow into the inner 
cavity. 

Reduction of rf, leads to a change of a constant M disk 
SED — its peak is now around 1 fim (thin solid line) and 
the peak of vF^ is about an order of magnitude higher 
than before. But the variation of the spectrum of a self- 
consistently evolved disk (thick solid line) is far more 
dramatic — it now peaks in the optical at 0.3 /im and 
the peak value of i^F^, is « 40 times higher than for a 
constant M disk. Between 1 /im and 500 fj,m the shape 
of the SED is well fit by vF^ oc z^^^/^, as expected for a 
constant Fj disk (SC95), now occupying the inner third 
of the radial extent of the disk and accounting for most 
of its luminosity. 

To summarize, the SED of a disk affected by the torque 
of a central binary is steeper, brighter, and extends to 
shorter wavelengths than the SED of its constant M 
counterpart having the same inner radius and mass ac- 
cretion rate Moo at large distances. 

These features (especially the steepness of the spec- 
trum) in principle make it possible to predict the exis- 
tence of a compact central binary in a disk based on the 
broadband spectroscopy alone, even in the absence of 
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Fig. 10. — Dependence of the disk spectrum on the fraction of 
mass inflow x = M{rin)/ Mao that penetrates across the orbit of 
the secondary and leaves the circumbinary disk. Different curves 
correspond to different x as labeled on the panel. Calculations are 
done for the situation depicted by solid curves in Figure|9l i.e. after 
4 X 10^ yr of evolution of the disk+binary system shown there. 



other indications of a SMBH binary such as the double- 
peaked line profiles caused by the relative motion of the 
binary components. This method may be the only way 
of inferring the presence of a binary in systems with face- 
on orientation. Relative brightness of disks affected by 
the binary torque should facilitate the detection of such 
systems out to large distances. 

5.3.1. Sensitivity of the SED to mass inflow across the 
secondary orbit. 

So far our calculations of the SMBH binary evolution 
and electromagnetic signatures have explicitly assumed 
that the binary torques completely prevent mass infiow 
from the circumbinary disk across the orbit of the sec- 
ondary. This allowed us to use the boundary condition 
M(rin) = in all our calculations, resulting in constant 
Fj disk near the binary. 

If the zero inflow requirement is relaxed and M(ri„) 7^ 
then one can still construct both the steady state and 
the evolving self-similar solutions as demonstrated in 
"i'Z 2.21 & 14 4 II Based on these solutions one can easily 
extend our results for the orbital evolution of the SMBH 
binary presented in !j5 5.1l fc l5 5.2l to the case of non-zero 
mass infiow at the inner edge of the circumbinary disk. 

Here we only show how the SED of the disk evolves as 
one varies the transparency of the inner barrier presented 
to the gas inflow by the binary torques. In Figure [TOl we 
show disk spectra computed assuming different values 
of X = M{rin)/Moo < 1- Instead of equation we 
now approximate the spatial distribution of Fj by the 
following simple formula: 



Fj{r) = 



Moo [I 



infl 



- X ihnti 



r > nnfl, 

- l{r))] , r < n„fl, 



(68) 



where l{r) = (GA^cr)^/^ ^infl = (GMcn^fi)^^^ . This pre- 
scription consists of two steady state solutions continu- 
ously matched at r = r-mti but with different AI inside 
and outside of this point. For x = this formula natu- 
rally reduces to equation ([M)) . 

Our calculations assume a system for which a spectrum 
is shown in Figure[H]at t = 4x 10^ yr for M(ri„) — 0, i.e. 
an equal mass, Mc = 10^ Mq SMBH binary at rf, = 10"^ 
pc surrounded by a disk with Moo = 10~^MEdd, and 
Hnfi ~ 0.03 pc. We show the spectrum of only the cir- 
cumbinary disk, i.e. in this calculation we do not account 
for the emission produced by an accretion disk(s) around 
the primary and/or secondary which should form when 
mass flows across the orbit of the secondary. 

One can see that as the transparency of the tidal bar- 
rier X increases towards unity the disk spectrum steadily 
approaches that of a constant M — Moo disk. This is not 
at all surprising because in the limit x ^ 1 disk structure 
reduces to that of a constant M disk, see Figure O As a 
result the slope of the power law portion of the spectrum 
steadily goes down from 12/7 to 4/3 as x is varied from 
to 1. 

It is also clear from Figure [10] that SED is strongly 
affected compared to the case of a constant M (x = 1) 
disk even if only a small amount of inflowing mass gets 
accumulated at the inner edge of the disk. For example, 
peak amplitude of the disk spectrum for x = 0.98 (im- 
plying that only 2% of the accreting mass gets stopped 
by the binary torques) is a factor of 2 higher than in the 
X = 1 case. And if the tidal barrier allows penetration 
of only 30% of the gas across the gap, the SED of the 
circumbinary disk is hardly distinguishable from that of 
a X = disk with no gas infiow at ri„. 

These results imply that the broadband SED of the 
disk is a rather sensitive measure of even a small amount 
of matter penetrating into the cavity cleared by the 
SMBH binary. Coupled with the measurements of the 
SED of the accretion disk(s), which may form around 
each of the binary components if x 7^ 0, these observa- 
tions can inform us on the efficiency of the binary torques 
at clearing a clean cavity at the center of the system. 

6. SUMMARY. 

In this work we explored the coupled evolution of a 
SMBH binary and a gaseous disk around it. Disk prop- 
erties (surface density, temperature, etc.) evolve under 
the action of binary torques, which constrain the flow 
in the inner part of the disk. To study this problem we 
have re-formulated evolution equations in terms of the 
angular momentum flux Fj. This significantly simplifies 
treatment of the steady state disk structure, when Fj is a 
simple linear function of the specific angular momentum 
I. 

We derived the disk properties as a function of Fj in 
different physical regimes that may be realized in cir- 
cumbinary disks around SMBH binaries. We demon- 
strated that radiation pressure can limit the value of Fj 
in disks around massive SMBH binaries by making the 
disk geometrically thick and susceptible to launching a 
radiation-driven wind. 

When the external mass supply to the disk at large 
distances is not matched at the inner edge of the disk be- 
cause of the binary torques, the disk evolves towards es- 
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tablishing a quasi-steady state in the inner region, where 
the local viscous timescale is shorter than the evolution 
time of the system. Viscous angular momentum flux in 
the inner disk steadily grows in time, which accelerates 
orbital evolution of the binary. 

We explored the dependence of this general picture on 
the system parameters (mass of the binary, mass accre- 
tion rate through the disk, etc.) and found the following 
in agreement with previous studies. 

1 . Tidal coupling to a circumbinary disk can substan- 
tially (by orders of magnitude) shorten the lifetime 
of the binary (IPP; Lodato et al. 2009; Haiman et 
al. 2009). 

2. For a long period before the GW emission takes 
over, the binary evolves in the limit when the mass 
of the secondary is much larger than the local disk 
mass (Haiman et al. 2009). 

3. Disk-driven evolution of the binary can be mea- 
surable by the space based gravitational wave an- 
tennae for low q systems with relatively low Ale 
(Kocsis etal 2011; Yunes et al. 2011). 

4. Spectrum of the disk affected by the binary torques 
is different from that of a conventional constant M 
disk: it extends to shorter wavelengths and more 
power is emitted. Instead of j^F^ cx ly^^^ the SED of 
a circumbinary disk exhibits a power law segment 
with vF^ cx (SC95). 

We also obtain a number of new results, summarized 
below. 

1. Self-consistent evolution of the disk resulting in a 
pile-up of mass at its inner edge accelerates the 
orbital evolution of the binary. 

2. Disk-binary coupling has a non-local character: the 
torque acting on the binary is determined by the 
state of the disk far from the binary, at the radius 
of influence nnfl, which steadily increases in time. 



rather than by the disk properties in the immediate 
vicinity of the binary. 

3. Evolution of the binary orbits exhibits a phe- 
nomenon of hysteresis — dependence of the evo- 
lution on the past history of the system, which is 
caused by the non-locality of the disk-binary cou- 
pling. 

4. Radiation pressure can strongly affect the disk 
structure even in cases when the mass accretion 
rate at large distances (in the constant M portion 
of the disk) is considerably sub-Eddington. 

5. Gas overflow across the orbit of the secondary 
affects binary mainly (or only) during the GW- 
dominated phase of its orbital evolution and is most 
important for low q, high Moo systems. 

6. Spectra of circumbinary disks strongly depend on 
the ability of accreting gas to cross the orbit of 
the secondary or otherwise leave the system, thus 
giving rise to a non-zero value of M at the inner 
edge of the disk. 

This list clearly implies that properly accounting for 
the fully self-consistent, time-dependent evolution of cir- 
cumbinary disks is crucial for understanding gas-assisted 
SMBH mergers. This general conclusion will hopefully 
inspire re-evaluation of some of the existing results for 
the orbital evolution of SMBH binaries and their ob- 
servational manifestations, both in the electromagnetic 
and the GW domains. Results of this work can also be 
extended to studying circumbinary disks around stellar 
mass binaries. 

I am indebted to Pavel Ivanov, Bence Kocsis and 
Zoltan Haiman for careful reading of the manuscript, 
open exchange of opinions, and a number of useful sug- 
gestions. The financial support for this work is provided 
by the Sloan Foundation, NASA grant NNX08AH87G, 
and NSF grant AST-0908269. 
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APPENDIX 

SCALING RELATIONS FOR ARBITRARY POWER LAW OPACITY. 

Here we summarize scaling relations for disk properties that result when opacity is a power law function of gas 
temperature T and density p (here taken to be represented by their midplane values): 



(Al) 



We make two additional assumptions regarding disk properties: (1) disk is optically thick and (2) radiation pressure 
is negligible compared to the gas pressure (radiation pressure dominated case is described by equations (P ^ - psp ). 
Combining equations ([T5)) . & (|26l) . and v = ac^/f2 valid in the gas pressure dominated regime one finds 



E(r): 
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where e = 10 + 3pi — 2^2- 
Using these results we can derive an expression for the diffusion coefficient Dj in the form (|5T 
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See Lyubarskij & Shakura (1987) for similar results. 



(A5) 
(A6) 



SUMMARY OF THE DIFFUSION COEFFICIENT BEHAVIOR IN DIFFERENT REGIMES. 



Here we summarize the behavior of the diffusion coefficient Dj in the power law form (I5ip and of the self-similar 
exponent n defined by equation (j56p for different objects and in different regimes explored in this work (see also 
Lyubarskij & Shakura 1987). 

In the radiation pressure dominated case f ^3 3.1 3.1.ip one finds for & = 



and for 5=1 
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Expression (jB2l) also holds true for the gas pressure dominated case with k — k^s f ^3 3.1 3.1.2p . 
In the gas pressure dominated case with k = Kff (SJ33II3II31) 



2-5 kq ^k 
3tt^ a \fi 



15/2 



a 



1/10 



3 

10' 



(B3) 



